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ABSTRACT 

\ 

The  objective  of  this  grant  is  to  pursue  research  in  the  field  of  noise  driven 
nonlinear,  dynamical  systems  by  means  of  both  analog  and  digital  simulations.  i  Topics 
for  these  simulations  are  normally  chosen  from  contemporary  theoretical  or /experi¬ 
mental  works,  usually  in  consultation  with  the  appropriate  groups.  The^eifiphasis  of 
the  research  is  on  laser  applications,  but  other  topics,  often  more  generally  related  to 
laser  dynamics,  have  also  been  undertaken.  In  the  area  of  laser  dynamics,  two  pro¬ 
jects  have  been  begun  and  one  is  still  in  progress>»The  first,  is  an  electronic  circuit 
modelling  of  the  effects  of  purely  multiplicative  noise  on  the  correlated  spontaneous 
emission  laser,  carried  out  in  consultation  with  W.  Schleich  (Max-Planck,  Munich), 
H..  Risken  andTUs-greup  (University  of  Ulm)  and  M.  O.  Scully  (Max-Planck  and 
University  of  New  Mexico)^^The  second.'''^ich  is  still  in  progress,^  is  a  modelling  of 
the  general  problem  of  stochastic  resonance,  and  a  study  of  the  effects  of  the  modu¬ 
lation  and  noise  as.  they  relate  to  the  recently  observed  phenomenon  in  a  dye  ring 
laser.  In  related  arejisT^four  projects  have  been  completed;  The  first  is  a  generaliza¬ 
tion  of  mean-first-passage^  time  calculations  to  problems  involving  spatio-temporal 
noise.  Specifically  we'Tiav^considere(^a  Brownian  particle  moving  in  a  random  spa¬ 
tial  potential  driven  by  temporal  noise,  and  we  highlight  an  unsolved  problem  which 
appears  to  be  extremely  difficult,  hut  "nevertheless  is  of  wide  physical  applicability, 
i  Second,  wo-hctv^  complete^  an  experiment  and  a  contemporary  theory,  oh  generalized 
switching  processes  in  the  presence  of  colored  noise. 'Though  the  object  of  this  work 
was  to  study  the  noise  behavior  of  a  generic  switchj there  is  an  obvious  application 
to  switching  in  lasers  and  optically  bistable  system^  This  work  was  carried  out  in 
collaboration  with  C.  Van  den  Broeck  (Limburg^  University  and  the  University  of 
Brussels)  and  P.  Hanggi  (University  of  Augsburg).4>Third.  -we  -have  -made"  a  discovery 
using  analog  techniques  with  strongly  colored  noise:  a  noise  correlation  time  induced 
change  in  the  topology  of  the  two  dimensional  probability  density  of  multistable  sys¬ 
tems.  ^e  unmistakable  implication  of  this  observation  for  first  passage  time  and 
bifurcaticnr  theory-is  that  (he  traditional  methods  which  reduce  the  problem  to  an 
"equivalent"  one  dimensional  form  are  grossly  inadequate,  and  that  the  multidimen¬ 
sional  nature  of  such  problems  must  be  fully  accounted  fo'r>Fourth.  we  have  com¬ 
pleted-  an  interesting  example  problem,  on  state  dependent  diffusion  posed  by  Lan- 
dauer  and  Van  Kampen.  While  not  offering  any  really  new  knowledge,  it  highlights 
an  old  question  which  has  frequently  been  ignored  especially  in  the  context  of  evolu¬ 
tion  theories.  Finally,  a  considerable  effort  has  been  expended  in  the  preparation  of  a 
review  article  on  colored  noise  dynamics  which  has  been  solicited  by  D.  ter  Haar  for 
Physics  Reports.  We  anticipate  that  this  article  will  be  finished  in  the  Spring  or 
early  summer  of  1989.  The  coauthors  are:  P.  Hanggi.  P.  Jung,  and  H.  Risken.  Con¬ 
currently,  new  data  acquisition  and  analysis  equipment  has  been  purchased  in  order 
to  modernize  the  simulation  laboratory;  a  process  which  is  not  yet  complete. 
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INTRODUCTION 


In  this  report,  progress  under  the  grant  named  on  the  cover  sheet  is  summarized 
for  the  period  12/01/87  to  11/31/88.  Six  projects  were  completed  in  the  sense  that 
these  were  works  in  six  identifiably  separate  topical  areas  and  that  each  has  lead  to 
a  publication.  The  topical  areas  spanned  a  range  from  first  passage  time  problems  to 
stochastic  resonance  and  state  dependent  diffusion.  They  are  described  in  more  detail 
below  and  in  the  attached  papers  in  the  Appendix. 

Each  of  these  projects  represented  a  distinct  application  of  the  technique  of 
analog  simulation.  Specifically,  that  means  that  an  electronic  circuit  model  of  a  par¬ 
ticular  system  was  constructed,  and  driven  by  a  noise  voltage  of  well  known  statisti¬ 
cal  properties.  Measurements  of  the  response  of  the  circuit  under  various  conditions 
were  made  and  compared  to  contemporary  stochastic  theory.  Specific,  physically 
observable  quantities  of  interest  to  the  theorists  were  measured:  the  first  passage  time 
and  its  probability  density;  stationary  one-  and  two-dimensional  probability  densities; 
time  evolving  probability  densities;  and  power  spectra  and  correlation  functions.  Since 
we  are  dealing  with  stochastic  systems,  each  of  these  quantities  must  be  appropriately 
averaged. 

The  experimental  process  is  as  follows:  An  electronic  circuit  model  of  the  system 
of  interest  is  designed  and  built.  The  "systems"  themselves  were  invariably  repre¬ 
sented  by  one  or  a  set  of  Langevin  equations.  Noise  from  a  noise  generator  is  passed 
through  a  linear  filter  to  establish  its  correlation  time  and  then  applied  to  the  circuit 
model.  The  fluctuating  responses,  i.e.  the  noisy  output  voltages,  were  digitized  and 
then  processed  by  a  computer  in  order  to  average  the  quantities  of  interest.  These 
measurements  are  then  analyzed,  the  error  bars  established,  and  then  compared  to 
theory  in  the  same  way  that  experimental  data  is  treated.  In  this  way,  analog  meas¬ 
urements  frequently  point  the  way  for  experiments  on  the  actual  physical  systems 
and  can  often  identify  any  pitfalls  which  might  arise  in  the  course  of  the  measure¬ 
ments.  This  latter  characteristic  has  been  especially  useful  in  the  experiments  on  sto¬ 
chastic  resonance  where  instrumental  band  width  played  a  crucial  role  and  was  not 
immediately  recognized  in  either  the  experiments  or  in  previous  digital  simulations.* 

The  remainder  of  this  report  is  organized  as  follows:  the  description  of  the  six 
completed  projects  is  followed  by  a  list  of  the  personnel  who  have  worked  on  the 
projects  and  the  degrees  granted  based  on  research  work  completed  under  the  grant. 
Certain  University  contributions  to  this  research  as  well  as  other  grant  support  is 
outlined.  Then  follows  an  outline  of  projects  currently  in  progress.  The  special  prob¬ 
lems  in  which  we  have  an  interest  and  the  objectives  of  these  current  projects  are 
discussed.  After  this  we  outline  our  plans  for  the  second  year  of  the  grant.  Two 
areas  here  seem  to  be  particularly  attractive  and  timely.  As  is  our  usual  practice, 
discussions  with  the  appropriate  theorists  have  taken  place  or  are  planned.  It  is 
important  to  emphasize,  though,  that  I  work  in  a  rapidly  developing  field  where  cer¬ 
tain  "targets  of  opportunity"  sometimes  rapidly  appear.  The  projects  on  stochastic 
resonance  and  noise  correlation  time  induced  topologies  were  just  such  examples.  It 
is,  therefore,  not  prudent  to  propose  a  rigid  structure  for  our  research.  The  last  sec¬ 
tion  contains  an  outline  of  our  budget  expenditures  for  both  personnel  and  equipment 
and  our  projections  of  these  for  next  year.  Finally,  reprints  and  preprints  of  papers 
which  were  produced  with  support  of  this  grant  are  attached  in  the  Appendix. 
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SUMMARY  OF  RESEARCH  PROJECTS  COMPLETED 

1)  First  passage  time  with  spatio-temporal  noise.  The  mean  first-passage  time 
has  for  many  years  been  a  very  useful  concept  for  studies  on  a  wide  variety  of 
problems  from  chemical  reaction  rates  to  theories  of  evolutionary  biology.  Often  inter¬ 
esting  applications  are  found  in  various  types  of  disordered  systems,  yet  in  only  one 
previous  instance  has  an  attempt  been  made  to  evaluate  the  MFPT  in  a  realistic 
random  spatially  distributed  potential.^  By  realistic,  we  mean  a  non-discontinuous 
random  functional,  as  would  be  found  in  macroscopic  natural  systems,  in  contrast  to 
the  usual  distribution  of  delta  functions**^.  Here  we  consider  a  classical  Brownian 
particle  moving  in  a  one-dimensional  random  potential  U(x)  which  has  a  non  zero 
correlation  length  f,  that  is.  colored  spatial  noise.  The  dynamics  are  given  by  the 
Langevin  equation, 

X  -  -dU(x)/dx  +  |(t),  (1) 

where  the  correlations  of  the  random  functions  are  given  by 

<U(x)U(x’)>  »  (C/l)exp-lx  -  x'l/f,  (2) 

and 

<{(t)Kt’)>  -  2D5(t  -  V).  (3) 

Note  that  the  spatial  noise  of  intensity  C  given  by  Eq.  (2)  has  a  correlation  length  I, 
while  the  temporal  noise  of  intensity  D  given  by  Eq.  (3)  is  white. 

In  spite  of  intense  effort  by  the  theorists:  A.  Engel.  L.  Schimansky-Geier  and  P. 
Hanggi,  ail  of  whom  worked  on  this  problem  to  some  extent,  the  more  realistic  prob¬ 
lem  -  colored  spatial  and  colored  temporal  noise  -  could  not  be  effectively  treated  in 
any  known  approximation.  We  highlight  this  as  an  open  problem  in  the  published 
version®.  The  white  temporal  noise  problem  which  we  did  pursue  was  resolved  into 
four  expansions,  framed  in  terms  which  are  easy  to  interpret  and  apply  by  experi¬ 
mentalists.  which  apply  in  various  regimes  of  parameter  space:  small  or  large  tempo¬ 
ral  noise  intensity  compared  to  the  spatial  noise  intensity  in  the  two  length  regimes: 
small  or  large  separation  between  the  exit  boundaries  compared  to  the  spatial  correla¬ 
tion  length.  The  Einstenian  relation  for  Brownian  motion,  <x^  -►  Dt,  is  recovered  in 
all  appropriate  limits.  A  discussion  of  the  possibilities  for  observing  the  predicted 
higher  order  departures  from  the  Einstein  relation  is  included.  The  details  can  be 
found  in  the  reprint  in  the  Appendix. 

Though  considerable  effort  was  expended  on  analog  simulations  of  this  problem, 
no  reliable  results  were  obtained,  so  only  the  theory  paper  was  published.  TTie  reason 
is  that  we  were  unable  to  obtain  in  a  reasonable  amount  of  experimental  time  a  suit¬ 
able  ensemble  average  over  sufficiently  many  realizations  of  the  spatial  potential.  We 
have  ideas  for  advancing  this  technology,  i.e.  to  significantly  speed  the  spatial  aver¬ 
age  over  realizations,  which,  however  require  the  new  data  analysis  systems  just  now 
being  installed.  It  should  be  noted,  however,  that  this  work  resulted  in  an  extremely 
practical  advance  in  analog  simulation  technology:  the  ability  to  simulate  any  arbi¬ 
trary  spatial  potential  in  both  one  and  two  dimensions.  This  means  that  we  are  no 
longer  limited  to  simulations  of  Langevin  equations  with  polynomials  made  up  of 
powers  or  trigonometric  functions  of  x  and  y,  but  can  simulate  truly  arbitrary  func¬ 
tions  so  long  as  they  are  finite. 
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2)  Noise  correlation  time  induced  topologies.  This  project  was  completed  in 
collaboration  with  F.  Marchesoni  of  the  University  of  Perugia,  and  it  represents  an 
actual  discovery  which  subsequently  prompted  a  degree  of  theoretical  activity.  The 
observations  were  made  initially  on  a  special  "soft”  bistable  potential  which  has  been 
used  in  soliton-soliton  interactions: 

U(x)  -  -ln{chx/(ch*x  +  sh^x)}.  (4) 

which  increases  only  linearly  for  sufficiently  large  x.  [It  should  be  noted  that  this 
potential  is  illustrative  of  our  present  ability  to  simulate  an  arbitrary  function  as  out¬ 
lined  in  the  previous  section.]  We  simulated  the  system, 

X  -  -dU/dx  +  5(t),  (5) 

with  strongly  colored  noise, 

<|(t)^(t’)>  -  (D/T)exp-lx  -  x’1/t.  (6) 

"Strongly  colored",  means  that  the  dimensionless  correlation  time  t,  which  measures 
how  fast  the  noise  fluctuates  on  the  time  scale  of  the  electronic  model,  is  »  1  (it 
was  typically  3  to  7  depending  on  the  depth  of  the  wells  and  the  strength  of  the 
noise  D).  We  measured  the  joint  probability  density  P(x,i).  Previous  studies  using  the 
standard  quartic  bistable  potential,  U  -(l/2)x*  +  (l/4)x^,  had  shown  characteristic 
colored  noise  topological  features®,  notable  the  so-called  "skewing  effect"  which  were 
simultaneously  confirmed  by  Risken’s  group  using  matrix  continued  fraction  tech¬ 
niques^.  Though  the  shape  of  the  cross  sections  through  these  densities  continuously 
evolved  with  increasing  t,  the  topology  remained  the  same.  Specifically,  the  two 
maxima  (corresponding  to  the  two  potential  wells)  were  always  separated  by  a  single 
saddle  point  which  was  always  located  on  the  x  axis  (at  the  origin  for  the  standard 
quartic).  The  measurements  on  the  soft  potential,  however,  showed  a  striking  new 
feature:  at  a  certain  critical  value  of  t.  the  single  saddle  point  "exploded"  into  two 
symmetrically  oriented  but  skewed  saddles  located  well  off  the  axis®. 

Confirmation  of  this  phenomenon  was  soon  obtained  for  an  entirely  different 
potential  (in  fact  a  periodic  potential)  by  the  Risken  group.  And  later  the  same  phe¬ 
nomenon  was  observed  in  analog  simulations  of  the  standard  quartic  potential  and 
even  of  a  random  potential®.  So  why  had  this  feature  not  been  previously  observed 
in  either  digital  or  analog  simulations  or  with  the  matrix  continued  fractions?  The 
reason  is  that  a  shallow  well  depth  (as  was  the  case  for  the  "soft"  potential)  is 
required  to  observe  the  effect  for  moderate  values  of  both  D  and  t. 

The  implications  of  this  discovery  for  mean  first  passage  time  theory  are  signifi¬ 
cant.  Extant  theoretical  approximations  for  colored  noise,  the  so-called  "small  r"  theo¬ 
ries,  had  always  reduced  the  necessarily  higher  dimensional  Fokker-Planck  equation 
to  an  "equivalent"  one-dimensional  form.*®  This  amounts  to  choosing  a  most  probable 
path  for  the  motion  of  a  particle  from  one  well  to  the  other  which  is  a  straight  line 
crossing  the  saddle  point  at  the  origin.  The  measured  topologies  for  t  greater  than  the 
critical  value  demonstrated,  in  contrast,  that  the  most  probable  path  which  now  must 
pass  through  one  of  the  off-axis  saddles  was  a  wide  curve  and  sometimes  a  tortuous 
path  which,  however,  always  avoided  the  origin.  Obviously,  in  order  to  treat  first 
passage  time  problems  for  strongly  colored  noise  a  multidimensional  approach  is  nec¬ 
essary.  This  has  recently  been  developed  by  Matkowsky**.  In  any  case,  the  observa¬ 
tions  present  an  interesting  problem  in  bifurcation  theory  which  has  stimulated  the 
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interest  of  several  groups'^. 

3)  Generic  switching  in  the  presence  of  noise.  Noise  necessarily  affects  all 
macroscopic  (and  probably  most  microscopic)  switches.  As  computer  memory  element 
densities  tend  toward  ever  larger  numbers,  with  a  concomitant  reduction  in  the  size 
of  the  individual  switch  elements  toward  dimensions  limited  by  quantum  noise,  the 
subject  of  noise  immunity  becomes  increasingly  important.  Moreover,  switches  operat¬ 
ing  even  at  the  quantum  level  are  better  described  with  colored  noise  theories*^. 

In  this  project  and  the  following  one,  we  have  undertaken  to  study;  first,  the 
effects  of  colored  noise  on  a  classical  generic  switch;  and  second,  the  predicted  effect 
of  noise  color  on  a  quantum  system,  the  correlated  spontaneous-emission  laser  (CEL). 
In  both  cases,  it  has  only  recently  become  possible  to  handle  the  problem  of  colored 
noise  driven  switching  with  adequate  approximate  theories**-’®. 

In  this,  the  classical,  project,  we  have  considered  a  generic  switch  as  repre¬ 
sented  by  a  general  Langevin  equation  of  the  form, 

X  -  x[-l  +  A(t)/(1  +  x^)]  +  m  (7) 

driven  by  the  usual  colored  noise 

<4(t)^(t’)>  -  (D/r)exp-|t  -  t’l/T.  (8) 

TTie  forcing  as  shown  in  Eq.  (7)  may  seem  strange,  but  it  is  the  result  of  the  simplest 
"soft"  potential  which  is  still  nonlinear  and  bistable  could  be  written: 

U(x)  -  (l/2)[x®  -  A(t)ln(l  +  x^).  (9) 

[A  soft  potential  is  one  for  which  the  restoring  force  increases  only  linearly  as  x  -► 
oo.]  Switches  with  "harder"  potentials  can  easily  be  represented,  but  the  generic 
switch  should  be  represented  by  a  soft  potential  for  which  the  effects  of  the  noise 
will  be  greatest.  Moreover,  this  potential  is  a  one-dimensional  representation  for  the 
dynamics  of  a  ring  laser.  In  the  above  expressions,  A(t)  represents  the  switch  signal: 

A(t)  -  Ao[-l  +  2u(t)],  (10) 

where  u(t)  is  the  Heavyside  (or  unit  step)  function.  Whence,  A(t)  -  -Aq  for  all  t  <  0 
and  A(t)  -  +Ao  for  all  t  >-  0.  What  this  really  means  is  that  the  potential  is 
switched  (in  a  time  short  compared  to  all  other  time  scales  in  the  simulation)  at  t  -  0 
from  a  potential  which  is  globally  stable  at  x  -  0  to  one  which  is  bistable.  The  col¬ 
ored  noise  |(t)  is  always  present,  and  forms  a  random  initial  condition  which  drives 
the  switch  off  the  unstable  state  (on  top  of  the  barrier  separating  the  stable  states  for 
t  >  0). 

The  probability  density  and  the  first  moment  of  the  time  necessary  for  the 
switch  to  respond  by  crossing  a  predetermined  threshold  were  measured  and  calcu¬ 
lated.  This  is  a  generalization  of  the  more  frequently  encountered  (but  less  physically 
realistic)  problem  in  the  decay  of  an  unstable  state  wherein  the  initial  condition  is 
taken  as  a  delta  function  and  some  mean  first  passage  time  theory  is  applied.  The 
general  approximate  theoretical  framework  for  this  problem  has  been  developed  by 
Dhara  and  Menon**,  and  we  have  applied  it  specifically  to  the  generic  switch’®.  For 
more  details,  see  the  paper  reproduced  in  the  Appendix’®. 
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More  recently,  Sancho  and  San  Miguel  and  their  coworkers,  stimulated  by  this 
work,  have  reworked  and  expanded  on  the  same  problem  with  an  improved 
theory*®’'^.  We  regard  this  body  of  work  as  a  significant  addition  to  the  literature  on 
the  subject  of  the  decay  of  unstable  and  meta  stable  states  or  more  specifically:  noisy 
switching  processes. 

4)  Noise  quenching  in  simulations  of  the  correlated  spontaneous-emission 
laser  (CEL).  This  project  is  a  natural  outgrowth  of  the  previous  one.  Though  we  do 
not  consider  switching  in  this  simulation,  we  do  test  some  predictions  regarding  noise 
quenching  in  the  CEL  both  via  analog  modelling  and  by  matrix  continued  fraction 
techniques.  Colored  noise  theory  is  only  just  beginning  to  be  applied  to  quantum  sys¬ 
tems*®.  The  colored  noise  theory  seems  to  be  difficult,  since  for  this  particular  model 
the  noise  is  both  colored  and  multiplicative.  The  dynamics  of  the  quantum  system 
can.  however,  be  reduced  to  a  description  by  a  classical  Langevin  equation*®: 

Vif  -  a  -  bsinti  +  ^(t)sin(i///2)  (II) 

with  the  usual  colored  noise  given  by  Eq.  (8).  In  this  case,  it  was  a  discovery  of  the 
analog  simulation  that  noise  color  (that  is  to  say.  increasing  noise  correlation  time) 
further  enhances  the  noise  quenching  already  resulting  from  squeezing  in  the  CEL. 
This  result  has  since  been  corroborated  by  matrix  continued  fraction  methods  and  by 
a  quantum  colored  noise  theory*®,  and  may  have  significant  practical  implications  for 
the  design  of  reduced  quantum  noise  measuring  devices®**.  For  further  details,  see 
preprint  in  the  Appendix. 

5)  Stochastic  resonance.  This  project  is  an  example  of  how  rapidly  analog  sim¬ 
ulation  can  respond  to  arising  topics  of  current  interest.  Stochastic  resonance  (SR)  is 
the  name  given  to  the  phenomenon  of  noise  assisted  transitions  in  bi-  or  multi-stable 
systems.  It  is  a  purely  nonlinear  effect  whereby  the  signal-to-noise  ratio  of  a  switch-- 
ing  signal  can  be  considerably  enhanced  by  the  deliberate  addition  of  noise.  It 
fits  very  well  indeed  into  our  program  along  with  the  two  previous  projects  on  the 
study  of  noisy  switching  systems.  SR  was  advanced  about  a  decade  ago  as  a  possible 
explanation  of  the  observed  periodicity  in  the  recurrences  of  the  Earth’s  ice  ages®*-®®. 
An  interesting  experiment  wherein  SR  was  observed  in  a  dye  ring  laser  has  recently 
refocused  attention  on  the  problem®®.  I  only  learned  of  the  laser  experiment  from  R. 
Roy  at  the  Los  Alamos  workshop®  in  March  1988.  well  before  his  Phys.  Rev.  Letter®® 
was  written.  It  was  nevertheless  possible  to  complete  an  analog  simulation  of  this 
system  by  early  summer  and  to  prepare  a  publication*.  While  the  laser  system  is 
multidimensional  (at  least  two  dimensional),  the  basic  phenomenology  of  SR  can  be 
represented  by  just  two  Langevin  equations:  one  which  represents  additive  modula¬ 
tion. 


X  -  X  -  X®  +  fcoscjt  +  $(t)  (12) 

and  the  other  representing  multiplicative  modulation, 

X  -  [1  +  €COSWt]*/®X  -  X®  +  |(t).  (13) 

where  the  coswt  term  represents  the  periodic  modulation  and  $(t)  is  the  noise.  The 
forcing,  x  -  x®,  is  just  that  due  to  the  "standard"  quartic  potential.  U(x)  -  -(l/2)x®  + 
(l/4)x*.  In  the  first  case,  both  the  modulation  and  the  noise  are  additive.  TTiis  means 
that  the  potential  is  "tipped"  up  and  down  by  the  modulation  so  that  the  potential 
wells  alternately  are  raised  and  lowered  in  relation  to  the  barrier  separating  them.  In 


the  second  case,  the  wells  remain  stationary  while  the  height  of  the  barrier  is  modu¬ 
lated. 

In  both  cases,  the  probability  that  a  transition  from  one  well  to  the  other  will 
occur  is  also  a  periodic  function  of  time.  Thus  the  problem  is  non  stationary,  and 
hence  raises  interesting  theoretical  challenges'**^.  What  had  escaped  previous  notice 
in  both  digital  simulations  and  in  experiments  was  the  fact  that  the  power  spectrum 
of  this  system  (at  least  that  of  Eq.  (12))  must  contain  a  delta  function  at  the  modula¬ 
tion  frequency.  That  observation  is  crucial  for  all  experiments-  and  simulations.  It 
means  that  efforts  to  measure  the  signal-to-noise  ratio  by  measuring  the  amplitude  of. 
the  peak  in  the  power  spectrum  can  at  best  be  only  qualitative.  This  observation 
prompted  our  comment^  In  our  simulation,  the  "input"  and  "output"  signal-to-noise 
ratios  were  measured  for  the  two  systems  represented  by  Eqs.  (12)  and  (13).  Our  res¬ 
ults  were  compared  with  the  McNamara  and  Wiesenfeld  theory^®.  For  more  details 
see  preprint  in  the  Appendix. 

6)  State  dependent  diffusion.  This  simulation  did  not  contribute  any  new 
knowledge  but  rather  was  undertaken  at  the  request  of  Rolf  Landauer  in  order  to 
illustrate  and  emphasi2e  an  old  point:  that  the  location  of  the  point  of  global  stability 
in  a  multi  stable  system  can  be  changed  with  multiplicative  (or  "state  dependent") 
noise.  We  chose  an  interesting  example  for  which  three  independent  theoretical 
approaches  had  been  taken.®®"®*  For  this  example,  a  standard  quartic  potential  was 
chosen,  with  the  noise  intensity  everywhere  constant  except  at  one  discontinuous  "hot 
spot"  near  the  potential  maximum  but  located  on  one  side  of  the  barrier.  It  was 
something  of  a  challenge  to  simulate  this  "state  dependent"  noise  with  electronic  cir¬ 
cuits.  The  instantaneous  location  of  the  trajectory  was  sensed  with  comparators  and 
sample-and-hold  circuits,  while  some  "extra"  noise  was  switched  in  when  the  trajec¬ 
tory  was  located  within  the  "hot  spot".  We  mad  measurements  of  the  stationary  pro¬ 
bability  densities  as  a  function  of  the  magnitude  of  the  extra  noise  in  the  hot  spot. 
Our  results  were  compared  to  the  theories  of  Landauer  and  of  Van  Kampen  (except 
for  a  differing  definition  of  particle  velocity,  the  results  are  equivalent,  though  the 
approaches  are  quite  different:  Van  Kampen’s  being  totally  thermodynamical  while 
Landauer’s  was  totally  statistical.)  Together  these  three  papers®®-®'*®®  make  an  interest¬ 
ing  pedagogical  collection  on  the  particular  point  regarding  global  stability. 

Buttiker’s  paper®®,  on  the  other  hand,  has  more  practical  consequences.  He  pred¬ 
icts  a  "noise  induced  probability  current"  in  a  continuous,  spatially  periodic  one¬ 
dimensional  system  with  periodically  placed  hot  spots.  Thus  a  mechanical  current  can 
be  the  result  of  spatially  periodic  heating.  There  may  be  applications  of  this  theory 
to  the  technology  of  superlattice  fabrication.  We  have  also  simulated  this  system,  and 
observed  and  measured  Buttiker’s  probability  currents,  but  have  not.  yet  written  up 
those  results  for  publication.®® 


PERSONNEL.  CONSULTANTS  AND  DEGREES  GRANTED 


1)  Students.  Two  graduate  students  were  supported  during  one  summer  and  par¬ 
tially  supported  during  one  academic  year.  They  both  received  MS  degrees  in  August 
1988.  Their  names  and  thesis  topics  are: 

Mr.  Goutam  Debnath,  "Observation  of  stochastic  resonance  in  a  bistable  potential  as  a 
model  for  a  bidirectional  ring  laser."  MS  thesis;  University  of  Missouri  at  St.  Louis; 
July  1988. 

Miss  Kathakali  Sinha,  "State  dependent  diffusion  (SDD)  in  a  bistable  potential".  MS 
thesis;  University  of  Missouri  at  St.  Louis;  July  1988. 

Unfortunately,  these  two  students  did  not  elect  to  remain  at  this  University  to  pursue 
PhD  work. 

In  addition,  an  undergraduate  student  Mr.  Gabor  Schmera  was  supported  part 
time  during  the  summer  and  academic  year. 

2)  Visitor.  Professor  Ting  Zhou  of  the  Institute  for  Semiconductors.  The  Chinese 
Academy  of  Sciences,  Beijing  has  been  a  visitor  from  May  1988  to  May  1989.  He  has 
been  supported  completely  by  funds  made  available  from  the  University  of  Missouri. 
He  has  worked  full  time  on  the  projects  described  herein. 

3)  Consultants.  My  grant  provided  for  theoretical  consultants.  The  following  con¬ 
sultants,  listed  along  with  the  topics  they  were  consulted  on,  were  supported  to  some 
extent  by  the  grant: 

a)  Prof.  H.  Risken,  University  of  Ulm,  ERG 
Noise  correlation  time  induced  topologies. 

b)  Prof.  P.  Mandel.  Free  University  of  Brussels.  Belgium 
Noise  driven  bifurcations  in  lasers.  Chaotic  responses  of 
coupled  COj  lasers. 

c)  Prof.  P.  Hanggi,  University  of  Augsburg,  FRG 
Stochastic  resonance.  Review  of  colored  noise  in  double 
wells:  an  invited  review  article  in  preparation  for 
Physics  Reports. 

d)  Dr.  W.  Schleich,  Max-Planck  Institute  for  Quantum  Optics,  Munich 
Noise  quenching  in  correlated  spontaneous-emission  lasers. 

e)  Dr.  C.  Van  den  Broeck.  Limbergs  Universitair  Centrum  and 
University  of  Texas  at  Austin. 

Random  potentials.  Decay  of  an  unstable  state. 

f)  Prof.  J.  M.  Sancho,  University  of  Barcelona,  Spain 
Decay  of  an  unstable  state  driven  by  colored  noise. 
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UNIVERSITY  CONTRIBUTIONS  TO  THIS  RESEARCH  PROJECT 

Before  this  grant  was  awarded  the  University  of  Missouri  committed  matching 
support  to  a  maximum  of  $80,000.  To  date  only  $13,000  of  this  has  been  made  avail¬ 
able.  That  money  has  gone  entirely  for  the  support  of  Professor  Zhou  as  mentioned 
above.  I  have  been  assured,  however  that  the  full  amount  promised  will  be  made 
available. 
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PROJECTS  CURRENTLY  IN  PROGRESS 

1)  Noise  correlation  time  induced  topologies.  This  project  has  resulted  in  one 
Letter  publication®  and  one  poster®,  but  is  far  from  complete.  In  collaboration  with 
the  Risken  group  at  Ulm,  a  complete  study  of  this  new  critical  feature  is  being  car¬ 
ried  out  over  all  parameter  ranges.  We  expect  this  work  to  be  finished  in  the  spring 
or  early  summer  of  1989. 

2)  Stochastic  resonance.  This  project  has  resulted  in  one  Comment  which  is  also 
letter  length^,  but  is  also  not  complete.  In  addition  to  signal-to-noise  ratios,  that  is 
power  spectra,  we  hope  to  measure  a  different  quantity  which  must  also  be  sensitive 
to  the  correlation  between  the  output  and-  the  modulation,  the  density  of  residence 
times.  TTiis  quantity  has  never  been  measured  in  an  SR  system  nor  is  there  yet  a 
theory  for  it.  It  is  a  quantity  which  is  sensitive  to  time  correlated  changes  in  the 
transition  probability.  In  addition,  this  quantity  does  not  suffer  from  the  defect  of  the 
power  spectra  in  that  there  should  be  no  delta  functions  in  the  output.  We  also  wish 
to  measure  the  total  power  in  the  band  of  both  the  signal  and  the  noise  in  order  to 
test  a  remarkable  prediction  of  the  McNamara- Wiesenf eld  theory®®:  that  the  the  signal 
actually  absorbs  power  from  the  noise  near  the  maximum  in  the  signal-to-noise  ratio 
curve.  It  has  to  date  been  impossible  to  test  this  prediction  in  the  laser  experiment. 

3)  Review  article  for  Physics  Reports.  As  stated  above,  this  is  a  review  to  be 
coauthored  by  Moss.  Jung.  Hanggi  and  Risken.  We  expect  it  to  be  finished  by  early 
or  mid  summer  1989. 

4)  Time  evolution  of  probability  densities  in  oscillators  with  weak  nonline¬ 
arities.  We  are  just  at  the  beginning  of  this  project.  The  appropriate  Langevin  equa¬ 
tion  is 


X  +  7X  +  kx  -  ex®  -  0,  (14) 

where  7  is  the  damping,  k  the  spring  constant  and  e  is  a  (weak)  nonlinearity  parame¬ 
ter.  In  the  case  that  the  damping  is  noisy,  7  7o  +  |(t),  the  linear  (k  -  0)  version  of 
this  oscillator  is  known  to  explode  (all  moments  -♦  00)  for  some  value  of  the  noise  in¬ 
tensity.  It  is  not  known  how  weak  nonlinearity  can  affect  this  behavior.  Moreover, 
for  the  first  time  ever,  the  time  evolution  of  the  probability  density  P(x,t)  has  been 
obtained  exactly  for  the  linear  case.  This  result  is  due  to  J.  Wm.  Turner  and  R. 
Lefever  of  the  University  of  Brussels  and  has  not  yet  been  written  up  for  publica¬ 
tion.  We  have  already  constructed  a  simulator  of  this  oscillator  and  begun  the  meas¬ 
urements  of  P(x.t).  The  oscillator  with  noisy  damping  is  itself  not  physically  interest¬ 
ing,  but  the  explosion  can  be  viewed  as  a  universal  bifurcation  event.  (The  usual 
normal  form  analysis  in  conventional  bifurcation  theory  always  begins  with  a  lineari¬ 
zation,  valid  for  small  parameters,  near  the  bifurcation  point.)  The  recent  solution  of 
the  time  dependent  Fokker-Planck  equation  for  this  linear  system,  may  lead  to  a  truly 
nonlinear  bifurcation  theory. 

5)  Postponements  and  advancements  of  noisy  Hopf  bifurcations.  This  is  a 
topic  which  we  have  been  interested  in  for  some  time.  We  have  previously  done 
experiments  on  both  th».  stationary®®  and  time  dependent®®  properties  of  the  Hopf 
bifurcation  in  the  Brusselator.  Of  interest  is  whether  noise  can  produce  an  advance¬ 
ment  or  a  postponement  of  the  Hopf  bifurcation.  The  original  analog  simulation  in 
the  Brusselator  showed  only  a  postponement®®.  One  theory  predicts  only  postpone- 
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ments^®.  while  two  others  predict  advancements^'^.  In  spite  of  our  earlier  measure¬ 
ments,  we  are  gathering  evidence,  from  precision  digital  simulations,  that  there  are 
indeed  advancements.  It  is  possible  that  the  advancements,  if  present  in  the  earlier 
analog  simulations  were  not  visible  due  to  a  lack  of  precision  combined  with  the  fact 
that  the  Brusselator  limit  cycle  is  highly  irregular.  The  digital  simulations  and  the 
theory  are  currently  being  applied  to  the  Hopf  bifurcation  in  the  Poincare’  oscillator: 

X  +  [A  +  B(x*  +  x^)]x  +  7X  +  ^(t)  -  0  (15) 

where  the  noise  can  either  be  additive,  as  shown  is  the  above  equation,  or  multipli¬ 
cative,  in  which  case  B  -*  Bq  +  |(t).  This  oscillator  has  the  advantage  that  the  limit 
cycle  is  regular,  i.e.  it  is  a  circle,  and  the  magnitude  of  the  velocity  of  the  trajectory 
x(x)  s  const.  This  question  is  important  in  the  analysis  of  real,  macroscopic,  physical 
systems  which  are  always  noisy.  Are  bifurcations  advanced  or  postponed  by  this 
noise?  And  if  so.  how  does  the  magnitude  of  the  advancement  (postponement)  depend 
on  the  ratio  of  the  noise  correlation  time  and  the  limit  cycle  period?  We  do  not  have 
good  answers  to  these  questions  yet  either  from  the  experimental  or  theoretical  points 
of  view  or  from  simulations.  Digital  data  are  slow  in  coming.  Near  the  bifurcation 
point  very  large  amounts  of  CPU  time  (the  order  of  hours)  on  our  VAX  8600  is 
required  just  to  obtain  one  two-dimensional  probability  density  P(x.x).  The_  analog 
simulations  are  much  faster  (order  of  minutes)  but  of  less  precision. 


PLANS  FOR  THE  SECOND  YEAR 

Obviously,  the  first  priority  is  to  finish  the  current  projects  listed  above.  How¬ 
ever.  there  are  two  additional  projects  which  we  would  like  to  pursue  time  and 
supply  of  students  permitting. 

1)  Weak  chaos:  weakly  dissipative,  weakly  nonlinear  oscillators.  This  project 
has  not  taken  form  yet  but  the  idea  is  to  examine  nearly  Hamiltonian  systems  using 
our  electronic  modelling  techniques.  We  are  following  here  the  work  of  G.  M.  Zas¬ 
lavsky  and  his  coworkers  regarding  this  possibility.  He  has  observed  the  initial  stages 
of  the  breakup  of  KAM  tori  at  the  onset  of  Hamiltonian  chaos  in  various  weakly 
driven  oscillators.  Using  digital  simulations,  we  have  been  able  to  reproduce  some  of 
his  results,  but  the  requirements  on  CPU  time  are  enormous.  Therefore,  if  it  is  possi¬ 
ble  to  study  any  of  the  phenomenology  using  analog  simulations  this  would  be  of 
great  interest.  In  particular,  2^1avsky  has  observed  the  formation  of  "stochastic 
webbs"  for  very  weak  chaos.  These  are  Just  chaotic  phase  plane  trajectories  close  to 
the  separatrices,  for  example  in  extended  systems  with  periodic  potentials  which  are 
driven  periodically^®"*^  Sometimes  the  webbs  display  interesting  symmetries,  e.  g. 
quasi  crystal  symmetries.  An  example  is  the  following: 

X  +  Wo^sinx  -  esin(kx  -  wt)  (16) 

For  rational  ratios  cj/cjo-  and  for  certain  (small)  values  of  c,  the  webb  is  formed 
around  the  phase  plane  trajectories  x(x).  Now  we  have  made  simulations  of  similar, 
though  non-Hamiltonian,  systems: 

X  +7X  +coo^sinx  •  ^(t)  (17) 


and  we  have  been  able  to  approach  Hamiltonian  behavior  for  7  0.  The  noise  ^(t), 

then  drives  the  system,  and  for  very  small  7  only  a  very  small  noise  intensity  is  nec- 
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essary  to  keep  the  system  on  a  trajectory  near  the  Hamiltonian  separatrix.  In  this 
case,  we  have  measured  the  two-dimensional  probability  density  P(x.x).  The  contours 
of  constant  probability  in  the  phase  plane  then  look  very  much  like  the  separatrix. 
We  call  these  densities  the  stochastic  phase  portraits*^.  In  fact  they  look  very  much 
like  Zaslavsky’s  stochastic  webbs.  but  of  course  our  system  is  not  chaotic  being 
driven  by  (weak)  high  dimensional  noise  instead  of  a  periodic  time  function.  It  cer¬ 
tainly  is  easy  to  build  an  analog  simulator  of  Eq  (16).  albeit  with  the  added  damping 
term  7X.  The  question  is,  by  careful  design  is  it  possible  to  make  7  small  enough  to 
unmistakably  observe  the  chaotic  webbs?  Alternatively  one  can  ask  if  any  of  the 
symmetries  observed  by  Zaslavsky  are  preserved  when  high  dimensional  noise  drives 
the  system?  The  actual  system  to  be  simulated  would  be  represented  by; 

X  +  7X  +a)o^sinx  -  €sin(kx  -  cjt)  +  |(t).  (18) 

It  is  probable  that  what  would  be  observed  would  be  some  mixture  of  the  chaotic 
webbs  and  the  noise  driven  portraits  with  5(t)  being  the  smallest  possible  inherent 
electronic  noise.  Could  these  yield  any  useful  information?  The  zaslavsky  group  have 
agreed  to  collaborate  with  me  on  these  simulations. 

[For  the  record  Michael,  we  contacted  the  Zaslavsky  group  nearly  8  months  ago, 
based  on  one  of  their  publications  in  Physics  Letters,  long  before  the  Physics  Today 
article  and  before  I  had  seen  the  Nature  article.  They  are  also  writing  a  book  for 
our  Cambridge  Nonlinear  Dynamics  series.] 

2)  Coupled  Oscillators:  a  model  for  two  coupled  lasers.  This  is  a  three  way 
loose  collaboration  between  myself  and  P.  Mandel  (Brussels)  who  is  doing  the  theory 
and  P.  Glorieux  (Lille)  who  is  doing  the  experiment  with  coupled  CO2  lasers.  They 
would  like  to  have  a  "catalogue  of  the  phenomenology"  from  us  in  order  to  point 
them  toward  regions  in  parameter  space  where  interesting  behavior  is  expected.  The 
system  is  five-dimensional,  and  is  expected  to  exhibit  complex  behavior.  We  have 
previously  modeled  a  single  laser  which  was  driven  by  an  external  oscillator  (instead 
of  a  second  laser)  and  observed  a  rich  dynamics,  including  a  Farey  Tree  of  phase 
locked  regions,  slow  mode  modulational  instabilities  (later  observed  in  real  lasers). 
Arnold  tongues  and  transitions  to  chaos”,  so  the  extension  to  two  coupled  lasers 
would  probably  be  straight  forward.  The  single  laser  model  is: 


E  -  -kE  +  (g"/7o)EN 

(19a) 

N  -  -7,(N  -  No)  -  (4g2/7o)|El2N 

(19b) 

with  periodically  modulated  losses 

k  -►  ko[l  +  mcos(wt)]. 

(20) 

where  E  is  the  field  intensity,  N  is  the  population  inversion,  the  7  are  the  two  polar¬ 
ization  modes,  g  is  the  atomic  field  coupling  factor,  and  k  is  the  cavity  loss.  Tne 
proposed  coupling  of  two  lasers  would  be  accomplished  by  coupling  the  cavity  losses, 
so  that  k,  -  bkj,  where  the  subscripts  now  represent  laser  1  and  laser  2.  The  system 
is  five-dimensional  with  two  variables  for  each  in  E,  N,  and  k. 

I  wish  to  add  a  caveat  here:  this  project  is  expected  to  be  long  and  difficult.  It 
will  necessitate  a  very  careful  cataloging  of  the  expected  wide  variety  of  complex 
dynamic  behavior.  I  cannot  complete  this  and  also  meet  my  other  commitments  with- 
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out  either  some  new  (and  very  good)  students  or  a  postdoctoral  student,  preferably 
both.  In  view  of  the  University’s  commitment  outlined  above,  money  is  not  a  prob¬ 
lem.  however  we  have  had  difficulty  in  attracting  good  graduate  students  here  rec¬ 
ently. 

3)Targets  of  opportunity.  It  is  necessary  to  save  both  some  time  and  laboratory 
capacity  in  order  to  respond  to  interesting  situations  which  may  arise.  During  this 
pjist  year  the  project  on  stochastic  resonance  was  just  such  an  opportunity. 
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The  mean  first-passage  time  of  a  particle  moving  in  a  spatially  correlated,  Gaussian  random  po¬ 
tential  and  subject  to  Gaussian  white  noise  is  calculated.  Depending  on  the  parameters  of  the 
problem,  we  find  four  different  regimes  of  diffusional  transport. 


The  mean  first-passage  time  (MFPT)  has  for  many 
years  been  a  most  useful  concept  in  the  analysis  of  rate 
processes  in  stochastic  systems.  While  the  origins  of 
MFPT  theory  are  traceable  to  early  Russian  work,'  more 
recent  pedagogical  treatments  are  to  be  found  in  Stratono- 
vich,^  Weiss,  ^  Van  Kampen,''  and  Risken,'  and  today  it 
continues  as  a  topic  for  active  research^  and  even  con¬ 
troversy.  ' '  In  these  previous  studies,  various  applications 
of  the  MFPT  theory  have  been  put  forth,  but  always  to 
bistable,  metastable,  or  periodic  potentials.  In  this  short 
note,  we  extend  the  MFPT  concept  to  include  particle 
motion  in  one-dimensional  random  fields  with  nonzero 
correlation  length  and  with  5-correlated  random  forcing; 
white  temporal  noise  driving  a  colored  spatial  noise  dy¬ 
namics.  We  hope  to  highlight  a  problem  in  contemporary 
stochastic  theory:  temporal  colored  noise  driving  a  spatial 
colored  noise  dynamics.  Though  such  a  problem  is  physi¬ 
cally  more  realistic  than  the  one  presented  here,  the 
theory  seems  to  be  more  difficult,  especially  since  all  tem¬ 
poral  colored  noise  theories,  even  for  completely  deter¬ 
ministic  potentials,  are  necessarily  approximate.*"" 

Stochastic  processes  in  random  fields  represent  an  in¬ 
teresting  problem  in  statistical  mechanics.  On  the  one 
hand,  there  is  a  variety  of  applications,  for  example,  in 
disordered  solid-state  systems  ^  or  models  of  biological 
evolution.'^  On  the  other  hand,  several  challenging 
mathematical  aspects  arise  in  their  theoretical  descrip¬ 
tion,'^"'*  and  interesting  applications  including  1//  noise 
arise. '*•*’  Previous  studies  mainly  concentrated  on  the 
stationary  probability  distribution  or  the  time  dependence 
of  the  moments  and  the  autocorrelation  function.  We  em¬ 
phasize  also  that  in  most  cases  previously  studied,  the  spa¬ 
tial  potential  was  5  correlated'^  in  contrast  to  the  nonzero 
correlated  one  used  here.  A  notable  exception  is  the  early 
work  of  dc  Gcnnes. 

We  consider  the  overdamped  motion  of  a  particle  in  a 
potential  Uix)  subject  to  additive  noise  ^(/).  Denoting 
the  position  of  the  particle  by  x(/),  the  equation  of  motion 
reads 

dtx{t)^—dxU{x)-¥^l.t),  (1) 

Given  two  points  x\  and  xi,  we  can  define  for  any  U{x) 
the  average  with  respect  to  ^(/)  of  the  time  at  which  the 
particle  for  the  first  time  reaches  x  2  if  it  was  at  xi  at  /  “O. 


This  is  the  well-known  MFPT  (Ref.  4),  T[[/(x)J,  which 
still  functionally  depends  on  L/(x).  Assuming  now  that 
U(x)  is  a  sample  of  a  random-field  ensemble  character¬ 
ized  by  some  functional  probability  density  Pll/(x)],  the 
question  about  the  ensemble-averaged  MFPT 

f-  fdU(x)PlL/(x)]TlU(x)]  (2) 

becomes  relevant. 

In  the  present  note,  we  consider  only  the  most  simple 
case  of  one-dimensional  motion  in  a  Gaussian  random 
field  with  Gaussian  white  noise,  i.e.,  we  take 

(^(r)>-0,  (5(/)^(/')>-2D5(/-/'), 

and 

P[C/(x)]  “■^exp  ( —  T  J*  j'dxdxV(x)A(x,x')L/(x')  , 

(3) 

where  D  measures  the  intensity  of  the  temporal  noise,  N  is 
a  normalization  constant,  and  A(x,x')  is  the  matrix  in¬ 
verse  of  the  spatial  correlation  function  B(x,x')  of  the 
random  field,  i.e., 

J'dx'A (x,x')B(x',x")  “5(x  — x") . 

Using  the  well-known  expression  for  the  MFPT  of  one¬ 
dimensional  processes  driven  by  Gaussian  white  noise,  ^  we 
find,  for  every  realization  f/(x), 

Tll/(x)]^-^fJdy-^fJdy'pHy') ,  (4) 

where 


p^Cxj—cxp 


t/(x) 


(5) 


denotes  the  stationary  solution  of  the  Fokker-Planck 
equation  corresponding  to  Eq.  (l).  Here,  and  in  the  fol¬ 
lowing,  we  assume  without  loss  of  generality  that  xi  “0 
and  X2  >  0. 

Equations  (4)  and  (5)  combine  to  give 


rl[/(x)}-jf^"dy'exp 


-jlU(y)-U(y')J 


(6) 
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Since  L'(.v)  occurs  in  (6)  only  linearly  in  the  exponential,  the  configurational  average  (2),  with  the  distribution  (3),  in¬ 
volves  just  a  Gaussian  functional  integral,  which  is  easily  performed  yielding 


For  macroscopically  homogeneous  random  fields,  we  have 
Bix.x')  =  5(1  .t  —  ,x'  I  ),  and  defining  5  =  5(0)  we  get 

^F'Jo  ^^v'exp  ^[5-5(>  -y’)]  ,  (8) 

or  equivalently. 

^•;^JoVX'^/v'exp  -^[5-5(>-')l  .  (9) 

Usually,  the  correlation  function  B(x)  is  assumed  to  fall 
off  exponentially  for  distances  large  compared  with  the 
correlation  length  /.  A  typical  choice  is 


5(.r)=-yexp(— |.x  I//) ,  (1 

1 -^(^  — l-F-^ 

2D  D-  X  2D^ 


where  C^Bil)  defines  the  intensity  of  the  static  disorder. 
In  the  following,  we  use  (10)  to  discuss  some  limiting 
cases  of  the  general  result  as  given  by  Eq.  (8). 

There  are  two  dimensionless  parameters  in  the  problem. 
The  first  one  is  x//  which  tells  us  whether  there  are  many 
variations  of  the  potential  U(x)  between  0  and  x  or  not. 
The  sedond  one  is  BlD^,  which  represents  the  ratio  of  the 
intensities  of  spatial  disorder  and  temporal  noise.  Four 
limiting  cases  are  of  interest  as  shown  in  Fig.  1. 

We  first  consider  the  case  5/D  '«1,  i.e.,  the  typical 
fluctuations  of  the  random  field  are  very  much  smaller 
than  the  intensity  of  the  temporal  noise.  Expanding  the 
exponential  in  (9),  we  find 


As  expected  (cf.  Fig.  1)  this  is  just  an  expansion  around  the  case  of  pure  diffusion,  which  is  characterized  by  the  familiar 
Einstein  result,^' 

(x-)-~Dt.  (12) 

Note  that  the  corrections  depend  on  the  ratio  x/l.  Forx<C/,  we  get  from  Eq.  (11) 

which  describes  the  influence  of  smooth  variations  in  the  random  field  on  T.  For  x  » /,  Eq.  ( U )  gives 

w  _  r  . 1  _  1  r  r  _  1 3  n 


=  _  X 

2D 


1+4  ^jL 

X  x^  ID* 


i, -31  +  14 

X  1x2 


+  0  e 

D^ 


which  accounts  for  the  corrections  of  T  due  to  many  but 
small  fluctuations  of  the  spatial  randomness  (see  Fig.  1). 

Of  more  interest  is  the  case  S/D^»a),  where  the  sp  <- 
tial  disorder  dominates.  For  x<K/,  we  then  have 
exp{(l/D^)l5-5(y/')l}=exp{5>'7D^/},  and  we  obtain 
from  (9), 

f=-^y-[exp(5x/D^/)  -  1  —  {BxIdH)]  . 

B  ^ 

For  very  large  values  of  5/D^  such  that  5x/D^/»  1,  we 
therefore  find  that 

f=(D7V5^)exp(5x/D^/).  (15) 

For  x«:/  the  average  variation  of  the  random  potential 
between  0  and  x  is  of  order  5x/ /,  and  the  exponential  fac¬ 
tor  in  Eq.  ( 1 5)  is  just  the  averaged  Arrhenius  factor  which 
controls  escapes  over  such  a  hill  in  the  presence  of  thermal 
noise  of  intensity  D. 

Finally,  we  consider  the  case  5/ D  1  and x'^l.  Now 
the  potential  shows  many  variations  between  0  and  x,  and 
the  particle  will  therefore  encounter  many  hills  and  val¬ 
leys  on  its  way.  First,  we  perform  the  y '  integral  in  Eq. 


0  1 

IXg-XJ  IL 


FIG.  1.  Schematic  plot  of  the  random  field  U(x)  between  two 
points  .ti  and  .t2  in  different  regions  of  the  parameter  plane. 
The  double  arrow  indicates  the  strength  of  the  temporal  noise. 
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(9)  exactly  and  obtain 


-£i 

B 

D- 

D- 

.(16) 


where  £i(x)  is  the  exponential  integral.  ■"  Because 
fil-x)-— oo  for  -v— *0,  the  main  contribution  to  the  in¬ 
tegral  in  Eq.  (16)  comes,  for  .x—  <»,  from  large  values  of 
y.  Hence,  we  can  neglect  the  constant  term,  and  using’’ 

£  I  (x )  ”  —  y  —  ln.v  —  X  —  ~  for  x  —  0 , 

n-i  nxnl 

(17) 

we  find  that 

f-(xV2£»)exp(£/Z)-)  .  (18) 

As  in  Eq.  (15),  the  main  contribution  comes  from  the  ex¬ 
ponential  factor,  which  has  a  simple  meaning.  On  its  way 
from  0  to  X  the  particle  has  to  climb  many  hills  (see  Fig. 
1),  and  the  height  (/max  of  each  is  a  Gaussian  random 
variable  w'ith  distribution 

£(£max)~exp(-L’L,/2S),  (19) 

as  follows  from  Eq.  (3).  The  time  r(t/max)  to  climb  such 
a  hill  is  given  by  the  Kramers  formula 

r((/max)-exp((/max/Z)).  (20) 


From  Eqs.  (19)  and  (20),  one  realizes  that  the  main  con¬ 
tribution  to  T  comes  from  fluctuations  with  height 


(21) 


Maxima  with  Umax^B/ D  are  quickly  climbed,  those  with 
Umax'^Bl D  are  very  unlikely  to  occur.  From  Eqs.  (20) 
and  (21),  we  get  T --exp(£/D^)  as  in  (18). 

The  X  dependence  in  Eq.  (18)  is  a  little  surprising  since 


the  result  corresponds  to  a  process  of  normal  diffusion 
with  renormalized  diffusion  constant  D—*Dexpi  —  B/ 
D^).  Usually,  one  expects  that  T  is  dominated  by  the 
transition  time  of  the  largest  fluctuation  between  0  and  X|. 
Using  the  distribution  Eq.  (19),  however,  one  realizes  that 
the  height  of  this  largest  fluctuation  varies_with  x  as 
(Inx)’^’  which  implies  a  contribution  to  T  of  order 
exp(Inx) For  large  x,  this  is  smaller  than  the  pure 
diffusion  time  T— x’ “exp(21nx)  and  Eq.  (18)  is  there¬ 
fore  explained. 

Note  that  the  problem  studied  in  Ref.  15  is  different 
from  ours,  insofar  as  there  the  derivative  dxUix)  was  tak¬ 
en  as  the  stochastic  variable.  The  corresponding  potential 
U(x)  was  therefore  a  Wiener  process  (also  having  zero 
correlation  length)  and  hence  it  has  statistical  properties 
very  different  from  ours  (see  Eq.  (3)1. 

In  conclusion,  we  have  studied  the  ensemble-averaged 
MFPT  for  the  overdamped  motion  of  a  particle  in  a  spa¬ 
tially  correlated  Gaussian  random  field  subject  to  Gauss¬ 
ian  white  noise.  The  general  results,  Eq.  (7)  and  (9),  have 
been  shown  to  interpolate  between  several  limiting  cases 
accessible  to  intuitive  arguments  and  possibly  to  experi¬ 
ment  as  well.  For  suflficiently  weak  spatial  disorder,  we 
found  expansions  around  the  normal  diffusive  behavior 
[Eqs.  (13)  and  (14)1,  i.e.,  the  MFPT  is  mainly  determined 
by  the  distance  to  be  traveled.  For  strong  spatial  disorder, 
the  order  of  magnitude  of  the  MFPT  is  given  by  an  Ar¬ 
rhenius  factor  corresponding  to  the  typical  fluctuation  of 
the  random  field  to  be  climbed  by  the  particle,  but  in  prac¬ 
tice,  in  this  limit,  the  MFPT  is  determined  by  the  local  de¬ 
tails  of  Uix). 

Discussions  with  Dr.  L.  Schimansky-Geier,  Professor 
W.  Ebeling,  Professor  P.  Hanggi,  Dr.  W.  Schleich,  and 
Dr.  F.  Marchesoni  are  gratefully  acknowledged.  This 
work  was  supported  in  part  by  the  U.S.  Office  of  Naval 
Research,  Grant  No.  N00014-88-K-0084. 
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Measurements  of  the  two-dimensional  probabihtv  density  using  an  analogue  simulator  of  a  system  w  ith  a  "soft"  bistable  poten¬ 
tial  show  a  new,  noise  induced  topological  feature:  the  appearance  at  large  noise  correlation  time  of  a  symmetric  pair  of  off-axis 
saddle  points.  The  implication  of  this  observation  for  mean  first  passage  time  calculations  is  discussed. 


The  problem  of  the  rate  of  escape  in  bistable  po¬ 
tentials  driven  by  colored  noise  has  recently  at¬ 
tracted  the  attention  of  a  number  of  investigators  ( for 
a  collection  of  reviews  see  ref.  [  1  ] ).  Let  us  consider 
a  one-dimensional  bistable  potential  l’(-v).  with  bar¬ 
rier  located  at  .v  =  0.  In  th  Smoluchowski  limit  the 
stochastic  rela.xation  of  t.ie  spatial  variable  is  de¬ 
scribed  by  the  following  equation  of  motion; 

.v=-r’(.v)-H((0  .  (1) 

The  simplest  choice  for  the  fluctuating  force  t{i)  is 
given  by  a  gaussian  noise  of  zero-mean  value  with 
autocorrelation  function  (acf) 

<t(Oe(0)>  =  {D/T)  exp(-U|/T)  .  (2) 

In  view  of  Doob’s  theorem  [2],  eqs.  ( 1 )  and  (2) 
can  be  replaced  by  a  markovian  system  of  two  sto¬ 
chastic  equations; 

.i=-V{.x)  +  e. 

i=(\/r)[-e(i)+i(n]  ,  (3) 

where  q{t)  denotes  a  zero-mean  valued,  gaussian. 
white  noise  with  acf  <i;(f)4(0)>  =2D<f(r).  Two 
properties  of  system  (3)  make  the  analytical  deter¬ 


mination  of  the  rate  of  escape  an  arduous  task;  ( i ) 
no  manifest  energy  detailed  balance  holds  [3j;  and 
(ii)  the  exact  stationary  distributions  p(.ir)  and  p(x. 
t )  are  unknown. 

The  very  definition  of  the  rate  of  escape  may  give 
rise  to  some  misunderstanding.  The  authors  of  ref 
[4].  for  instance,  identify  the  rate  of  escape  in  bi¬ 
stable  potentials  in  the  small  noise  limit  (D«  AT) 
with  the  smallest  non-vanishzng  eigenvalue  of 
the  corresponding  Fokker-Planck  equation.  This 
definition  is  particularly  well  suited  to  computations 
by  means  of  numerical  algorithms  based  on  contin¬ 
ued  fraction  expansions  [5].  In  the  case  of  analogue 
or  digital  simulations  A(t)  can  be  determined 
uniquely  as  the  reciprocal  of  the  slowest  decay-time 
of  the  temporal  acf  of  .v(f)  [4]. 

By  contrast,  much  effort  has  been  put  [6-10]  into 
characterizing  the  rate  of  escape  in  bistable  poten¬ 
tials  as  the  reciprocal  of  the  mean  first  passage  time 
[6]  (MFPT).  /"(r).  for  the  brownian  particle  to  cross 
over  the  potential  barrier.  In  the  white-noise  limit 
(t=0)  it  can  be  proven  that  the  r“'(0)  does  in¬ 
deed  coincide  with  A(0)  up  to  the  first  order  in  the 
small  parameter  D/W-  compare  the  results  of  refs. 
[11].  However,  the  equivalence  of  T  " '  (  t  )  and  A  ( r ) 
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is  suggested  not  to  hold  in  the  presence  of  colored 
noise  (t>0)  as  indicated  by  the  conflicting  predic¬ 
tions  for  the  T-dependence  of  the  rate  of  escape  in 
refs.  [4,7-10], 

MFPT  calculations  invariably  assume  that  for 
D<scAy.  the  escape  from  one  potential  well  into  the 
other  occurs  any  time  x  changes  sign.  In  the  two-di¬ 
mensional  space  (.r,  «)  this  means  that  the  barrier 
crossing  occurs  any  time  a  single  trajectory  intersects 
the  line  x=0.  the  most  probable  crossing  point  being 
located  at  (0,  0).  This  assumption  has  been  criti¬ 
cized  first  in  ref.  [10],  Subsequently,  the  authors  of 
ref.  [4]  determined  the  separairix  numerically,  that 
is,  the  curve  of  constant  probability  which  passes 
through  the  saddle  point  between  the  peaks  of  the 
distribution  function  p(.v,  e).  Since  p(x.  t)  cannot 
be  factored  into  a  simple  product  of  functions  of  t 
and  X.  it  is  nontnvial  to  determine  the  separatrix  for 
intermediate  to  large  values  of  r.  The  usual  answer 
to  such  a  criticism  is  that  in  the  limit  of  very  small 
values  of  T  [  12]  (but  how  small  is  hard  to  say)  the 
correlation  of  x  with  (  would  play  a  minor  role;  so 
that  the  most  probable  escape  path  crosses  the  po¬ 
tential  barrier  at  the  unstable  point  (0,  0)  as  in  the 
white-noise  limit.  The  effect  of  time-correlation  is 
then  supposed  only  to  modify  the  curvature  of  the 
effective  potential  at  the  bottom  of  the  well  and  the 
top  of  the  barrier.  Consequently,  the  corrections  to 
the  rate  of  escape  at  small  r  are  not  expected  to  de¬ 
pend  much  on  D.  at  least  for  small  D. 

In  this  Letter  we  present  a  result  of  an  analogue 
simulation  experiment,  which  illustrates  the  impor¬ 
tance  of  the  (X.  <)-space  structure  for  a  full  under¬ 
standing  of  the  relaxation  dynamics  of  the  class  of 
systems  described  by  eqs.  ( 1 )  and  (2).  We  chose  an 
unconventional  potential  function  Vix), 

which  is  motivated  by  the  analysis  of  the  double  sine- 
Gordon  kink  [13.14].  This  is  not  the  potential  used 
in  the  works  cited  in  refs.  [7-10],  however,  it  has 
the  great  advantage  that  certain  exact  results  can  be 
obtained  by  means  of  supersymmetric  transforma¬ 
tions  [14].  This  potential  is  shown  in  fig.  1  for  /? = 2. 
Numerical  values  of  the  relevant  potential  parame¬ 
ters  are;  the  barrier  height. 


2  - 


Fig.  1.  The  soft  potenlial  given  b>  eq.  (4). 


Ar=ln(:^^^)«0.67  . 

\2shRj 

the  positions  t.Xo  of  the  local  minima. 

.v„  =  (  ch-  '  sh  /?i0.98  . 

the  curvature  of  the  bottoms  of  the  wells. 
K"  (  ±.x„)  =  1  -  l/sh'/? 5:0.92;  and  finally,  the  cur¬ 
vature  at  the  barrier,  I  r"{0)|  =  1  —  2/ch-7?*0.86.  .A 
unique  feature  of  this  family  of  potentials  is  that  they 
are  linear  for  large  l.tUi.e,,  K(.y)~1x:|  forx->±oo. 
A  good  determination  of  A(0)  has  been  obtained  by 
means  of  a  supersymmetric  transformation  [14]  and 
the  identity  T  " '  ( 0 )  =  A  ( 0 )  has  been  established  for 
this  potential.  In  addition,  the  slow  divergence  of 
V(x)  for  large  |.y|  required  the  introduction  of  a 
modified  MFPT  formula  [15], 

We  have  constructed  an  analogue  simulator  of  eqs. 
( 1  )-(4)  using  an  already  well  established  technique 
for  simulating  arbitrary  potential  functions  [16], 
Measurements  of  the  two-dimensional  statistical 
density  p{x,  t)  for  several  values  of  r  are  shown  as 
contour  plots  in  fig.  2. 

The  main  features  detected  by  Risken  and  co¬ 
workers  [4.5  ]  and  observed  in  previous  simulations, 
both  for  the  quartic  double-well  potential  [17]  and 
a  sinusoidal  potential  [  1 6  ] ,  are  reproduced  here;  ( i ) 
the  skewing  of  p(jc.  t)  which  increases  with  r,  and 
(ii)  the  separatrix  which  does  not  coincide  with  the 
axis  at  x=0  [17], 

Most  notably,  however,  a  novel  feature  of  p{x.  t ) 
shows  up  for  values  of  r  which  are  large  compared 
with  1 1 '"  ( 0 )  I  ~  V"  ( .v,i )  -  0.9.  For  r  large  enough  ( the 
magnitude  increases  weakly  with  D )  two  symmetric 
saddlepoints  in  the  distribution  located  at  (  v.„  fj 
and  ( —x^.  —  t,)  become  observable  with  our  appa- 
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Fig.  2.  Measured  statistical  densities  pl.x.  <)  for  D=0A:  (a) 
r=0.1.orquasi-white  noise,  with  amplitudes  9"00.  3400  and  1100 
aiihc  peaks,  separainxand  lowest  contour  respectively,  (b)  r=1.5 
where  the  crater  has  disappeared.  The  peak,  separatnx  and  low¬ 
est  contour  amplitudes  are  54700,  470.  and  430  respectively,  (c) 
r=  5.0  showing  the  crater  ( shaded  area )  and  the  off-axis  pair  of 
saddle  points  (double  arrows ).  The  amplitudes  of  the  peaks,  sad¬ 
dles  and  lowest  contours  in  the  bottom  of  the  crater  are  98700. 
1 10  and  20  respectivelv  The  amplitudes  are  in  arbitrary  units. 


raius.  This  means  that  for  large  enough  r,  the  system 
presents  two  symmetric,  most  probable  escape  paths 
connecting  the  potential  wells,  which  make  a  detour 
around  the  single  crossing  point  located  at  (0,  0)  for 
T=0,  and  run  instead  across  the  saddle  points. 

It  is  not  yet  clear  whether  this  observed  noise  in¬ 
duced  change  in  the  topology  of  p(x,  ()  occurs  ^t  a 
specific  threshold  value  r,,,.  or  whether  the  two  sad¬ 
dle  points  instead  move  out  from  (0,  0)  continu¬ 
ously  for  T>0. 

An  example  is  shown  by  the  contour  plot  in  fig.  2c 
for  t=5.  where  the  shaded  area  surrounding  (0.  0) 
depicts  the  crater  in  p(x.  t).  The  two  saddle  points 
are  marked  by  the  double  arrows,  which  lie  on  the 
measured  contour  closest  to  the  separatrix,  and  are 
located  approximately  at  .v,~.Vo/2  and  For 

comparison,  densities  at  r=1.5.  where  the  saddle 
points  have  just  disappeared,  and  at  r=0. 1,  which  is 
close  to  white  noise  are  shown  in  figs.  2b  and  2a 
respectively. 

A  simple,  physical  explanation  for  this  observa¬ 
tion  is  the  following:  the  locations  of  the  wells  of  the 
potential  under  study  can  be  approximated  by 
|,v±.Xo|  (apart  from  an  additive  constant).  The 
presence  of  a  constant  external  force,  for  example,  of 
magnitude  ~e  would  correspond  to  a  momentary 
tilting  of  the  potential  due  to  the  addition  of  the  ex¬ 
tra  term  .re.  In  the  neighborhood  of  €  =  e,w  ±  1,  the 
effective  potential  K{r)+.«  becomes  unstable  and 
the  trajectory  falls  into  the  single  potential  well.  This 
mechanism  obviously  applies  only  when  r  is  large 
enough  to  allow  the  trajectory  to  relax  into  the  new 
potential  well  before  e  changes  again.  If  this  picture 
is  appropriate,  the  noise  correlation  time  induced 
change  of  topology  is  a  general  effect  which  should 
be  observable  for  an  arbitrary  pair  of  wells  separated 
by  a  barrier. 

It  is  clear  that  the  effective  activation  energy  cor¬ 
responding  to  the  newly  observed  crossing  points  (.Xs, 
ej  and  (  -  .r,,  -ej  is  different  (though  perhaps  not 
much  different  in  this  example)  from  the  .Arrhenius 
factors  modified  for  use  at  large  values  of  r,  for  ex¬ 
ample  by  the  ansatz  of  Hanggi  et  al.  as  defined  in  ref 
[  8  ] .  We  conclude  by  pointing  out  that  any  theory  for 
colored  noise  driven  bistable  potentials  should  be  able 
to  reprodu  e  the  complexity  of  the  stochastic  dy¬ 
namics  in  the  two  variable  space  (.v.  e). 


324 


Volume  131,  number  6 


PHYSICS  LETTERS  A 


29  August  1 988 


We  are  grateful  to  W,  Schleich  for  an  interesting 
discussion  and  to  Hans  Risken  for  a  recent  discus¬ 
sion  wherein  he  informed  us  that  the  off  axis  saddle 
points  have  ver>  recently  been  observed  at  large  cor¬ 
relation  times  for  systems  with  periodic  potentials 
using  the  matrix  continued  fraction  technique.  This 
w  ork  w  as  supported  in  part  by  the  L'.S.  Office  of  Na¬ 
val  Research,  grant  No.  N00014-88-K-0084. 

References 

(  1  ]  F  Moss  and  P  E.  McCliniock.  eds..  Noise  in  nonlinear 
dsnamical  ssstems,  V'ol.  1  (Cambridge.  1988).  to  be 
published. 

|2  )  W  Horsthemke  and  R.  Lefeser.  Noise  induced  transitions 
I  Springer.  Berlin.  1984  i  eh.  S. 

[  3  I  F  Marehesoni.  Ph\s.  Lett.  A  H)1  (1 9Sb  ill. 

( 4  J  P.  Jung  and  H.  Risken.  Z,  Ph\s.  B  61  l  1983  )  367; 

Th  Leiber.  F  Marehesoni  and  H.  Risken.  Phss.  Rev.  Lett. 
59  (  1987)  1381; 

F  .Marehesoni.  Phys.  Rev.  A  36  ( 1987)  4050. 

[5]  H.  Risken,  The  Fokker-Planck  equation  (Springer,  Berlin, 
1986): 

F  Marehesoni  and  P  Grigolini,  J.  Chem.  Phys.  78  ( 1983) 
628'', 

(6J  R.L.  Stratonovieh.  Topics  in  the  theory  of  random  noise. 
Vol,  1  (Gordon  and  Breach,  New  York.  1963)  p.  98. 

[■'IP,  Hanggi.  F  Marehesoni  and  P  Grigolini,  Z.  Phys,  B  5o 
(  1984  1  333. 

[8|P,  Hanggi,  T.J.  Mroekowski,  F.  Moss  and  P.V.E. 
McClintock.  Phys.  Rev.  A  32  i  1985)  695. 


[9|  .M  Sancho.  F  Sagues  and  M.  San  .Miguel.  Phys.  Rc\.  A  33 
(  19Sb  1  3304: 

R.F.  Fox.  Phys.  Rev  A  33  (  1986)  46'’,  34  (  1986|  525:  3 
I  1988)  91  1; 

J.  .Masolixer  K..  Lindeiiberg  and  B.J.  West.  Phys.  Rc'.  A  34 
(  1986  I  235.  35  I  1987  i  3086; 

C.R  Doering.  P.S  Hagan  and  C  D  1  exerniore.  Phys  Rex 
Lett.  59  (  1987)  2129: 

•M.M  Dygas.  B.J.  Matkowsky  and  Z  Sehuss.  Sl.A.M  J  Appl. 
Math.,  to  be  published. 

P.  Jung  and  P  Hanggi.  Phys  Rex  A  35  i  1987 )  4464.  and 
preprint; 

K.  Lindenberg.  B.J.  West  and  J  Masolixer.  in:  Noise  in 
nonlinear  dynamical  systems.  Vol.  1.  eds.  F.  .Moss  and  P  V  E. 
McCliniock  (Cambridge.  1988  )  ch.  4.  to  be  published 

{ 10)  P.  Hanggi  and  P  Talkntr.  Phys.  Rex  A32(I985)  1934 
[  1  I  I  R.S  Larson  and  M.D.  Kostin.  J  (  hem.  Phys  69(19^8) 
4821: 

(J  EdholmandO.  Leimar.  Phxsiea  A  9,s  i  |9'ui  313. 

[12]  JM  Sancho.  M  ban  Miguel  S  L  Kat/ aiui  . I  D.  (luniiKi. 
Phys.  Rex.  A  2o  (  1982  i  I  589 

I  13]  D  K.  Campbell.  M.  Peyrard  and  P  Sudano.  Physica  D  19 
(  1986)  165. 

P.  Sodano.  M.  El-Batanouny  and  C  R  Willis.  Phys  Rex  B 
34(1986)  4936. 

[  14 ]  M.  Bernstein  and  L.S.  Brown.  Phys.  Rev.  Lett.  52  l  1984 ) 
1933; 

R.  Dull.  A,  Kharc  and  L'.P.  Sukhaime.  Am.  J  Phxs  56 
(1988)163. 

1 15)  F  Marehesoni.  P,  Sodano  and  M.  Zanetti.  MIT.  CTP  No. 
I54U  1  1987), 

(16)  K.  Vogel.  H.  Risken.  W.  Schleich.  .M,  James.  F,  .Moss  and 
P  V.E,  MeClintoek.  Phys.  Re'.  A  35  i  1 98" )  46,3. 

[  I""]  F.  Moss  and  P.V.E,  McClintock,  Z.  Phys.  B6I  i  1985  )  381. 


325 


PHYSICAL  REVIEW  A 


VOLUME  38.  NUMBER  9 


NOVEMBER  I,  1988 


Switching  in  the  presence  of  colored  noise:  The  decay  of  an  unstable  state 

Mark  James* 

Department  of  Physics,  University  of  Lancaster,  Lancaster  LAI  4YB,  United  Kingdom 

Frank  Moss 

Department  of  Physics,  University  of  Missouri  at  St.  Louis,  St.  Louis,  Missouri  63121 

Peter  Hanggi 

Institute  of  .Mathematics,  University  of  Augsburg,  D-8900  Augsburg,  Federal  Republic  of  Germany 

Christian  Van  den  Broeck 

Limburgs  Universitair  Centrum,  B-3610  Diepenbeek,  Belgium 
(Received  14  March  1988;  revised  manuscript  received  29  June  1988) 

Switching  events  are  studied  by  means  of  a  parametrically  operated,  fast  transition  from  a 
monostable  to  a  bistable  potential  In  the  continuous  presence  of  colored  noise.  The  problem  is  thus 
the  decay  of  an  unstable  state  with  random  initial  conditions.  We  calculate,  using  contemporary 
colored-noise  theory,  and  measure  by  analog  simulation,  the  relaxation  time  to  cross  a  reference 
boundary,  and  we  contrast  this  with  the  strictly  defined  mean  first  passage  time. 


1.  INTRODUCTION 

In  this  paper  we  study  a  switching  process  under  the 
influence  of  noise  which  is  common  to  a  class  of  parame¬ 
trically  activated  bistable  systems.  The  system  exhibits  a 
single,  or  monomodal,  state  Xq  while  inactive,  but  devel¬ 
ops  a  bimodal  potential  upon  receipt  of  a  parametric 
switching  signal.  The  bimodal  potential  is  developed  in  a 
time  which  is  short  compared  to  the  characteristic 
dynamical  response  time  of  the  system.  The  initial  state 
Xq  thus  becomes  an  unstable  state  at  the  instant  when  the 
switching  signal  is  received.  This  state  is,  however,  per¬ 
turbed  by  noise  which  drives  the  decay.  After  a  relaxa¬ 
tion  time  (  r>,  which  we  analyze  and  measure  herein,  the 
system  settles  into  one  of  the  bistable  states. 

Such  generic  switching  processes  were  first  proposed 
by  Landauer'  as  possible  zero-  or  low-energy  switches. 
They  have  been  used  for  quite  some  time  as  examples  in 
discussions  on  the  minimum  energy  dissipation  necessary 
for  measurement  and  for  information  transmission  and 
computation.’  In  view  of  this  interest  and  of  the  funda¬ 
mental  importance  of  the  switching  process  itself,  it  is 
relevant  to  study  the  dynamics  of  switching  in  the  pres¬ 
ence  of  noise,  which  is  inescapable  in  macroscopic  sys¬ 
tems. 

The  first  theory  and  measurements  on  such  noisy 
switching  events  induced  by  linearly  swept  parameters 
and  their  remarkable  noise-averaging  properties  were 
carried  out  by  Kondepudi  et  al.^  Noisy  switching  events 
in  lasers  have  been  studied  experimentally  by  Zhu  et 
and  analytically  by  Broggi  et  al.*^^'  and  later  simu¬ 
lated  with  analog  circuits.^  However,  all  of  these  studies 
were  of  switching  events  induced  by  variations  of  the  pa¬ 
rameters  on  time  scales  comparable  to  the  characteristic 
dynamical  response  time  of  the  system.  By  contrast  in 
the  present  study  the  switch  parameter  is  operated  by  a 
step  function  at  time  f  =  0,  thus  preparing  an  initial  state 
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which  later  decays  into  the  induced  bistable  states. 

The  problem  of  the  decay  of  an  unstable  state  in  the. 
presence  of  white  noise  has  been  well  studied  by  a  variety 
of  techniques.*’’  Among  the  more  familiar  and  useful 
methods  for  describing  the  onset  of  macroscopic  order 
are  the  time  evolution  of  the  variance  the  mean 

first  passage  time  (MFPT),  and  the  onset  of  bimodality  in 
the  probability  density.  The  MFPT  in  recent  years  has 
become  a  widely  used  tool  for  analyzing  such  stochastic 
processes.* 

Non-numerical  colored  noise  theories  are,  however, 
necessarily  approximate.  This  arises  because  of  the  extra 
variable  necessary  to  describe  noise  with  nonzero  correla¬ 
tion  time.  The  Langevin  equation  and  its  analogous 
Fokker-Planck  (FP)  equation  are,  consequently,  at  least 
two  dimensional.  Since  the  latter  equation  can  be  solved 
exactly  only  in  one  dimension,’  recent  years  have 
witnessed  a  veritable  explosion  of  approximative 
schemes, all  of  which  seek  to  reduce  the  FP  equation 
to  an  “equivalent”  one-dimensional  form,  and  an  increas¬ 
ingly  vigorous  debate  regarding  their  accuracy  and  appli¬ 
cability.'^’”"'’  With  the  exception  of  Refs.  12  and  14  all 
of  these  are  adaptations  of,  or  improvements  on,  an  ap¬ 
proximation  originally  put  forth  by  Stratonovich^®’^' 
based  on  expansions  valid  for  small  correlation  times. 
How  accurate  they  are  for  a  given  correlation  time  de¬ 
pends  on  the  application  and  often  the  noise  intensity  as 
well.  Here  we  present  measurements  for  r  as  large  as 
t=5  for  which  both  one-  and  two-dimensional  stationary 
probability  densities  have  been  studied.’  '’’’’ 

In  this  work  the  nonsmall  noise  color  introduces 
difficulties  at  two  levels.  In  order  to  calculate  the  station¬ 
ary  mean-square  fluctuation  intensity  in  the  initial  state, 
we  will  use  a  simple  linearization  which  is  expected  to  be 
accurate  for  small  noise  intensities.  More  detailed 
theories  such  as  those  cited  above  could  be  used  instead, 
but  the  linear  theory  is  sufficiently  accurate  for  the 
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initial-state  calculations.  We  then  analyze  the  transient 
relaxation  toward  the  stable  state  using  colored-noise 
modified  Suzuki  scaling*  following  a  recent  theory.^^ 
These  results  are  then  compared  to  measurements  made 
on  an  analog  simulator  of  an  example  switch. 

We  find  that  the  relaxation  time  is  a  decreasing  loga¬ 
rithmic  function  of  the  noise  intensity,  results  that,  for 
white  noise,  were  anticipated  much  earlier. In¬ 
creasing  the  noise  correlation  time  increases  the  relaxa¬ 
tion  time  but  does  not  have  a  large  effect. 

This  paper  is  organized  as  follows.  In  Sec.  II  the 
theory  is  developed  following  Refs.  12  and  23.  In  Sec.  Ill 
the  simulator,  its  operation,  and  the  measurement  tech¬ 
niques  are  described.  The  results  are  presented  in  Sec.  IV 
and  compared  to  the  calculations.  Conclusions  and  a  dis¬ 
cussion  are  presented  in  Sec.  V. 

II.  THEORY 

First  we  consider  the  stationary  dynamics  in  a  single- 
well  potential,  i.e.,  before  the  control-parameter-induced 
switching  event.  The  potential  is  given  by 

f/(x  )  =  |[x^— /I  ln(  1 -f-x^)]  ,  (la) 

and  the  forcing  is 

f(x)=-x[-\->r  A(t)/{\+xh]  ,  (lb) 

where  .4  ( r )  is  the  control  parameter.  As  shown  in  Fig.  1, 
the  potential  is  monostable  for  A  <  \  and  bistable  for 
A  >1.  Moreover,  the  width  of  the  monostable  potential 
depends  upon  |  /4  | .  Our  model  is  determined  by  the 
Langevin  equation 

x—f{x)  +  ^(t)  ,  (2a) 

where  |;(  f )  is  an  exponentially  correlated,  Gaussian  noise 
source  with  zero  mean: 

<^(r)|:(5)>=(Z)/r)exp(-  I  r-5  1  /r)  ,  (2b) 

where  D  is  the  noise  intensity  and  r  is  the  noise  correla¬ 
tion  time.  The  system  is  prepared  in  the  monostable  ini¬ 
tial  state  by  setting  Ait)=  Aq<1,  for  which  the  deter¬ 
ministic  solution  Xo=0,  is  globally  stable.  After  having 
been  prepared  in  this  initial  state  for  a  time  sufficiently 
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FIG.  1.  The  potential,  showing  a  monostable  initial  state  for 
/4o=  — 1.38  and  the  switched  bistable  states  at  the  values  of  A 
indicated. 


FIG.  2.  The  circuit  diagram  of  the  electronic  simulator.  The 
multipliers  and  the  divider  are  standard  chips  available  from 
Analog  Devices,  Norwood,  Mass. 

long  to  establish  stationary  statistical  conditions,  the  sys¬ 
tem  is  switched  at  f =0  so  that  A(t)=  A  >1  for  all  r  > 0. 
The  control  parameter  is  thus  a  step  function. 

We  first  consider  the  initial  state  for  which  stationary 
conditions  are  assumed.  In  order  to  proceed  it  is  neces¬ 
sary  to  obtain  the  stationary  probability  density.  In  fact, 
we  will  only  need  the  second  moment  (x^)5,=  (xo(t)) 
of  the  stationary  initial  distribution.  Since  both  the  ex¬ 
periment  and  the  theory  to  be  used  later  are  restricted  to 
small  values  of  D,  linearized  theory  will  be  sufficient  for 
calculating  (xo(r)).  We  begin  with  the  following  FP 
equation,  valid  for  small  D ; 

Mo  I  +i)x-|:]P(x,g,n| 

^P(x,i,t) 

d§  T 


D  3^  . 

+  ^-^P(x,^,t)  .  (3) 

In  the  steady  state,  the  equations  for  the  moments  read 
(  I  /lol -l-l)<xJ(T))-<x^>=0  ,  (4) 

(  I /lol +l+T-')<x|>-<^^>=0  ,  (5) 

-<|O  +  ^=0,  (6) 


FIG.  3.  Example  trajectories  measured  for  /lo=  — 1-38  and 
4  =  2  for  2J  =0. 1  and  t=  1.0.  The  bistable  potential  is  switched 
at  f  =0.  The  vertical  scale  is  in  volts  with  the  potential  minima 
at  ±1.0  V.  The  horizontal  scale  is  in  ms. 
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FIG.  4.  Examples  of  the  probability  density  of  the  relaxation 
time  from  which  the  (T)  were  obtained  for  r=I.O,  A  =2.0, 
— 1.38,  and,  reading  from  left  to  right,  D,  =  11.7X  10"’ 
V’,  (r)  =  2.4  ms;  Z>,  =0.20x  10"’  V’,  (7’)=4.4  ms;  and 
Z),  =  2.56x  10"*’  V’,  {  r)  =  6.6.  The  horizontal  scale  is  in  ms. 

from  which  we  conclude  that 

<.x5(r))=D|(l  +  I  ^0  i  )[l+r(l  +  I  ^0  1  ])-'  +  0(D^)  . 

(7) 

This  describes  the  stationary  dynamics  before  the  switch¬ 
ing  event  takes  place. 

For  t  =0*,  the  system  is  switched  to  a  bistable  state  as 
shown  in  Fig.  1  for  .4  >0,  and  the  state  at  .t=0  is  ren¬ 
dered  locally  unstable.  Since  the  switching  event  takes 
place  in  a  time  very  short  compared  to  all  other  time 
scales  in  the  problem,  the  decay  is  driven  by  the  noise  dy¬ 
namics  described  by  Eqs.  (6)  and  (7)  toward  the  two  local¬ 
ly  stable  states  Xj  =  ±  ,4  —  1  created  at  f  =0. 

The  simulation,  described  in  Sec.  Ill,  measures  the 
residence  time  T,  for  a  random  walker  x{t)  necessary  to 
cross  a  reference  value  =x,/2  for  the  first  time  after 
the  switch  function  is  activated  at  t  =0.  The  mean  of  this 
quantity  (  D  is  a  measure  of  the  time  scale  on  which  the 
system  assumes  macroscopic  order  and  is  closely  related 
to  (but  not  identical  with)  the  MFPT  on  the  same  inter¬ 
val.  The  residence  time  is  like  a  weighted  MFPT  wherein 
the  distribution  of  initial  conditions  is  a  Gaussian  whose 
second  moment  is  given  by  Eq.  (7).  By  contrast,  our 
switch  is  more  physically  realistic  in  supposing  that  the 


FIG.  5.  <  r>  in  ms  vs  D  in  V’  for  Aq=  —  1.38  and  A  =2.0. 
The  values  of  7  are  0.1,  •;  0.5,  X :  1.0,  o ;  and  4.9,  0.  The  solid 
lines  are  fits  to  the  data  and  the  dashed  line  is  a  predicted  result 
for  7=0.5. 


D(V2) 


FIG.  6.  (  T>  in  ms  vs  D  in  V’  for  Ao=  — 1.38  and  A  =3.0. 
The  values  of  r  are  0.1,  •;  0.5,  X ;  1.0,0;  and  4.9,  0.  The  solid 
lines  are  fits  to  the  data. 

noise  is  continuously  present,  as  it  is  in  all  real  switches. 
In  fact  it  is  possible  that  the  instant  the  unstable  state  is 
born,  the  trajectory  xq(i=0)  could  already  exceed  the 
reference  boundry  Xg.  For  small  D,  however,  in  the 
range  where  both  the  theory  and  the  simulation  are  accu¬ 
rate,  such  events  are  extremely  rare. 

The  time  for  decay  from  an  unstable  state  has  been 
studied  by  many  authors.*’^^'^'  As  shown  first  by  Kubo 
et  the  relaxation  time  exhibits  a  logarithmic  depen¬ 
dence  on  the  noise  intensity.  Very  recently,  Suzuki’s  scal¬ 
ing  theory*'^'  has  been  generalized  to  include  colored 
noise,^^  with  the  result 

(r>  =  -(l/2a)ln(Cj(xg>-|-£>[a{l+aT)]-’l)  ,  (8) 

where  a~f'{x  =0'*')=  A  —1,  and  where  C  is  a  constant 
that  depends  only  on  the  parameters  of  the  deterministic 
nonlinear  flow  /lx ),  i.e.,  on  Aq  and  A,  but  not  on  D  or  t. 
The  theoretical  value  of  C  depends  also  on  the  detailed 
definition  of  the  residence  time,  for  example,  as  the  time 
for  the  second  moment  {x^){t)  to  relax  to  some  refer¬ 
ence  value  X/j ,  or  the  time  required  for  P(x,f )  to  become 
bimodal,  or  the  strictly  defined  MFPT,  etc.  But  certainly 
C  is  of  order  unity  so  that  InC  is  of  order  zero. 


FIG.  7.  ( D  in  ms  vs  D  in  V’  for  /lo=  —2.40  and  A  =4.0. 
The  values  of  7  are  0.1,  •;  0.5,  x  ;  1.0,  O;  and  5.0,  0.  The  solid 
lines  are  fits  to  the  data.  The  dot-dashed  and  dashed  lines  are 
theoretical  predictions  for  7=  1.0  and  0.1,  respectively. 
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TABLE  I.  Calculated  and  measured  parameters  for  =  —  1. 38.  A  =2,  and  Bg  =  —0.03. 


r 

_  th«of 
ffl 

m 

6"-“'(r)-t-Bo 

4.9 

1.15 

1.11 

0.80 

0.77 

0.77 

1.0 

1.15 

1.14 

0.23 

0.20 

0.20 

0.5 

1.15 

1.13 

0.08 

0.05 

0.10 

0.1 

1.15 

1.12 

-0.11 

-0.14 

0.02 

In  our  experiment,  the  noise  color  enters  Eq.  (8)  in  two 
places:  First,  there  is  a  colored-noise  effect  on  the  initial 
state  Pq{x,t)  which  leads  to  {xq{t))  as  given  by  Eq.  (7). 
Second,  the  dynamics  of  the  decay  is  dependent  on  r  as 
shown  explicitly  in  Eq.  (8).  Both  effects  move  (  T)  in  the 
same  direction,  i.e.,  increasing  r  results  in  increasing 

<r>. 

Observing  that  Inx  =ln(  lOHogj^x,  we  can  recast  these 
results  in  the  convenient  form 


<  r>  =  —  m  logioD  -|-6(t)-|-5o  , 

(9) 

where 

In  10 

m  = -  , 

liA-\) 

(10a) 

and 

Bo  =  -ni  log,oC~0  , 

(10b) 

so  that  6(t)  describes  the  combined  dependence  of 
on  the  noise  correlation  time: 

b{T)=-m  log,o|(  I  Ag  I  -l-l)“‘[l-l-r(  I  Aq  !  -1-1)]“' 

-h(  A-l)-'tl-t-T<^-l)]-'l  .  (11) 

These  results  can  be  accurate  only  in  the  range  of  small 
D.  Indeed,  we  used  linear  theory  to  calculate  (xq). 
Furthermore,  Suzuki  scaling  is  expected  to  be  accurate 
only  for  small  D,  but  here  the  quantitative  limit  is  not 
known. 

In  the  following  sections  we  describe  an  analog  simula¬ 
tion  of  Eqs.  (1)  and  (2)  for  the  range  0. 1  <r<5  and  for 
10“*<D  <  10“'.  Since  as  shown  by  Eq.  (2b), 

the  range  of  D  corresponds  to  a  range  of  noise  voltage 
y„=i  of  somewhat  less  than  three  orders  of  magnitude, 
or  about  10  mV  to  a  few  volts,  which  corresponds  to  the 
usable  dynamic  range  of  our  analog  simulators. 

III.  ANALOG  SIMULATION 

We  have  constructed  a  circuit  model  of  Eqs.  (1)  and  (2) 
using  by  now  quite  standard  techniques.^^  The  schematic 
diagram  of  this  simulator  is  shown  on  Fig.  2. 

The  noise  voltage  F„(f),  is  obtained  from  a  Wandell 


and  Goltermann^^  wide  band  (  >  KX)  kHz),  Gaussian 
noise  generator.  In  order  to  create  colored  noise,  ( t )  is 
passed  through  a  linear,  single  pole  filter  with  a  transfer 
function  H(o})  =  [l  where  co  is  the  radian 

frequency  and  the  noise  correlation  time.  The  charac¬ 
teristic  response  time  of  the  simulator  is  established  by 
the  integrator  time  constant  r,  ,  as  shown  on  Fig.  2.  The 
dimensionless  noise  correlation  time,  as  it  appears  in  the 
theory  and  in  Eq.  (2),  is  the  ratio  T=r„ /r,.  Whenr«l, 
the  simulator  perceives  the  noise  as  quasiwhite,  whereas 
T>  1  marks  the  colored-noise  regime.  With  V„{t)  =  ^{t), 
Eq.  (2b)  defines  the  noise  intensity  D  =  r{  for  t  =s. 
The  mean-square  noise  voltage  is  the  measured 

quantity  in  our  simulation  which  defines  D.  In  this  ex¬ 
periment  we  always  maintained  D  <0. 30,  and  the  corre¬ 
lation  time  varied  over  the  range  0. 1  <  r  <  5. 

In  the  absence  of  noise,  the  discrepancies  between  the 
deterministic  steady  states  measured  on  the  simulator  are 
smaller  than  ±2.5%  when  compared  to  the  steady-state 
(i=0)  solutions  of  Eqs.  (1)  and  (2a).  Measurements  in 
the  presence  of  noise  are,  of  course,  subject  to  statistical 
errors  which  can  be  reduced  by  increasing  the  number  of 
samples  in  a  given  average.  In  this  experiment  the  sta¬ 
tistical  errors  are  estimated  (from  the  repeatability)  to  be 
~±5%.  The  largest  error,  and  the  most  difficult  to 
quantify,  is  systematic  and  shows  up  in  the  quasiwhite 
noise  end  of  the  range  of  r.  This  results  from  the  limited 
dynamic  range  and  bandwidth  of  the  simulator.  For  r 
small,  F„(r)  is  large,  and  an  increasing  number  of  its 
large-amplitude  excursions  in  the  wings  of  the  Gaussian 
are  clipped  as  r  is  decreased.  At  r=0. 1,  this  results  in 
discrepancies  between  our  measured  probability  densities 
and  white-noise  solutions  of  the  FP  equation  which,  in 
places,  are  as  large  as  -^20%.  In  this  experiment,  relaxa¬ 
tion  times  are  measured  for  which  it  is  difficult  to  esti¬ 
mate  the  systematic  errors  at  the  small-r  end. 

In  operation,  a  rectangular  wave  Ait),  operating  be¬ 
tween  the  voltages  Aq  and  A  was  applied  to  the  divider 
as  shown  in  Fig.  2.  The  frequency  of  this  wave  was  ad¬ 
justed  so  that  in  the  state  Aq  sufficient  time  was  allowed 
for  the  initial  probability  density  of  x  to  become  station¬ 
ary  (several  hundred  times  r, ).  The  wave  switches  from 
Aq  to  A  at  t  =0  and  at  the  same  time  the  data-analysis 


_ TABLE  11.  Calculated  and  measured  parameters  for  Ao=  — 1.38,  A  =3,  and  go=0.07. _ 

_r _ _ b^^ir) _ ^ 

4.9  0.58  0.58  0.63  0.70  0.70 

I.O  0.58  0.58  0.31  0.38  0.39 

0.5  0.58  0.58  0.20  0.27  0.31 

0.1  0.58  0.57  0.07  0.14  0.28 
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TABLE  III.  Calculated  and  measured  parameters  for  —2.4,  A  =4.0,  and  Bo=Q.  16. 


r 

„  theor 
ftl 

^.hcor(^)^5o 

5.0 

0.38 

0.39 

0.55 

0.71 

0.71 

1.0 

0.38 

0.39 

0.32 

0.48 

0.47 

0.5 

0.38 

0.39 

0.24 

0.40 

0.39 

0.1 

0.38 

0.38 

0.12 

0.28 

0.35 

system,  a  Nicolet-Lab  80  computer  and  digitizer  connect¬ 
ed  to  the  circuit  at  x{t),  was  triggered.  A  time  series  of 
typically  2000K  digitized  points  was  then  obtained,  and 
the  time  at  which  the  trajectory  first  crossed  the  thresh¬ 
old  Xn=x^/2  {x^  —  ±v'a~1  are  the  deterministic 
steady  states)  was  measured  and  stored.  Ten  example  tra¬ 
jectories  are  shown  in  Fig.  3,  where  the  threshold  is 
marked  by  arrows,  and  an  example  crossing  time  at  T,  is 
shown.  After  a  large  number  of  such  measurements  (typ¬ 
ically  lO'*)  the  computer  tabulates  the  mean-relaxation 
time  ( D  and  its  density  P{T).  Three  examples  of  the 
densities  are  shown  in  Fig.  4  for  three  values  of  D.  As  ex¬ 
pected  (  T  >  increases  as  D  becomes  smaller. 

IV.  RESULTS 

The  results  of  our  simulation  are  summarized  on  three 
graphs  and  compared  to  the  theoretical  predictions  in 
three  tables.  Figure  5  shows  our  measured  values  of  <  T ) 
versus  D  for  four  values  of  r  as  indicated  by  the  different 
symbols  for  /4o  =  — 1.38  and  .4  =2.0  V.  Each  set  of  data 
were  matched  by  least -squares  fit  to  the  equation 

{T)  =  -m\oz^oD  +  b  ,  (12) 

and  the  values  and  were  extracted.  These 

were  compared  to  the  results  predicted  by  Eqs.  (9)-(ll). 
As  shown  by  Eq.  (9)  the  entire  colored-noise  eflTect  is 
represented  by  the  constant  hir).  It  is  evident,  however, 
that  the  non-color-dependent  constant  Bq  is  not  negligi¬ 
ble.  We  obtained  a  value  for  Bq  by  matching  the  experi¬ 
mental  and  theoretical  results  at  r=4.9  and  5  (where  the 
simulation  is  most  accurate)  and  then  compared  the  pre¬ 
dicted  and  observed  values  in  Table  I,  where  the  h‘*''°'^(T) 
and  m*'’*®'  are  calculated  directly  from  Eqs.  (10a)  and 
(11),  and  is  to  be  compared  with  h‘*'*‘’'’(r)-|-So.  The 
results  are  shown  on  Table  I.  The  largest  discrepancies 
are  for  r=0.5  and  0.1  as  expected.  Even  so,  on  the  loga¬ 
rithmic  scale  of  Fig.  5  the  difference  between  Eq.  (9)  using 
the  theoretical  values  and  the  data  are  small,  as  shown  by 
the  dashed  line  which  is  to  be  compared  to  the  crosses 
t=0.5).  The  solid  lines  through  the  data  for  t=0.  1  and 
t=4.9  are  example  plots  of  the  results  of  the  least- 
squares  fits  of  the  data  to  Eq.  (12). 

Figure  6  and  Table  II  show  the  results  for  /4o  =  — 1.38 
and  .4  =3.0.  It  is  evident  that  while  Bq  has  changed  con¬ 


siderably,  the  discrepancy  at  t=0.  1  is  decreased.  Figure 
7  and  Table  III  show  the  results  for  4o  =  — 2.4  and 
y4  =4.  On  Fig.  7  we  also  show  the  theoretical  result  for 
r=1.0  as  the  dot-dashed  line  to  be  compared  with  the 
open  circles.  The  dashed  line  is  the  prediction  for  r=0. 1 
and  shows  the  largest  disagreement  with  the  measure¬ 
ments  (closed  circles). 

V.  CONCLUSIONS  AND  DISCUSSION 

The  agreement  between  the  calculated  and  the  mea¬ 
sured  results  in  this  work  is  satisfactory  considering  the 
approximations  necessary  to  achieve  a  colored-noise 
theory.  We  emphasize  further  that  the  strictly  colored- 
noise  contributions  to  the  theory  have  two  sources:  the 
decay  theory  of  Dhara  and  Menon^^  and  the  ansatz.^^ 
The  latter,  which  agrees  with  the  linearization  result  in 
Eq.  (7),  can  be  used  to  calculate  the  second  moment  of 
the  initial,  stationary  (r<0)  density,  an  application  for 
which  it  is  known  to  be  relatively  accurate.  The  relative 
importance  of  each  of  these  contributions  depends  on  /Iq 
and  A  as  shown  by  the  two  terms  in  brackets  in  Eq.  (11). 
For  small  ^nd  large  A,  (xq)  dominates  the  decay 
process  and  hence  the  ansatz  is  more  important,  while  for 
large  Aq  and  small  A  the  reverse  is  true,  and  the  non¬ 
linear  decay  process  dominates.  We  conclude  by  pointing 
out  that  colored-noise-driven  decay  of  unstable  states 
should  find  applications  in  a  variety  of  switching 
scenarios  most  notably  in  laser  dynamics  and  nonlinear 
optical  bistability. 
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Comment  on  Stochastic  Resonance 


G.  Debnath*,  T.  Zhou+  and  Frank  Moss 
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University  of  Missouri  at  St.  Louis 
St.  Louis.  MO  63121 

ABSTRACT 

A  recent  and  interesting  experimental  paper  [B.  McNamara.  K.  Wiesenfeld  and 
R.  Roy.  Phys.  Rev.  Lett.  60.  2626  (1988)]  has  refocused  attention  on  the  problem  of 
stochastic  resonance  by  presenting  measurements  of  the  signal-to-noise  ratio  (SNR)  of 
a  noise  driven,  periodically  modulated  bistable  ring  laser.  We  point  out  that  the  theo¬ 
retical  SNR.  as  defined  in  this  and  a  previous  work,  is  always  infinity,  because  addi¬ 
tive  modulation  leads  to  a  delta  function  in  the  power  spectrum  of  the  output.  Quan¬ 
titative  information  on  stochastic  resonance  is  contained  in  the  strength  of  this  delta 
function  relative  to  the  noise  background.  We  qualitatively  reproduce  the  SNR  data 
with  an  analog  simulator  using  a  standard  quartic  bistable  potential.  In  this,  as  in 
previous  experiments  and  simulations,  a  peak  in  the  observed  power  spectrum  is  a 
reflection  of  the  delta  function,  but  the  amplitude  of  the  peak  is  rendered  finite  (and 
hence  measurable)  only  because  of  the  finite  resolution  of  the  measurement  system. 


PACS  numbers:  5.40.  +  j.  2.50.  +  s.  42.65.Pc 
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The  phenomena  of  stochastic  resonance  was  first  investigated  by  Benzi.  Sutera 
and  Vulpiani*  and  later  suggested  by  C.  Nicolis*  and  by  Benzi,  Parisi,  Sutera  and 
Vulpiani*  as  a  possible  explanation  of  the  observed  periodicity  in  the  recurrence  of 
the  ice  ages.  In  this  model,  a  pair  of  stable  climate  states  separated  by  a  barrier  is 
imagined.  This  bistable  system  is  driven  by  noise  resulting  from  random  fluctuations 
of  the  solar  constant.  Using  reasonable  climate  models,  it  was  demonstrated  that  such 
fluctuations  could  trigger  switching  events  between  the  two  stable  states  on  time 
scales'  which  are  in  approximate  agreement  with  the  period  of  the  observed  recur¬ 
rences  (10’  yrs)*-’.  The  switching  events  would,  however,  be  uncorrelated  random 
occurrences  in  time.  In  order  to  explain  the  observed  periodicity,  a  modulation  of  the 
height  of  the  barrier  or  of  the  alternate  depths  of  the  potential  wells  was  introduced 
into  the  model.  The  resulting  switching  events,  while  still  randomly  occurring,  must 
now  be  correlated  with  the  periodic  forcing  since  the  switching  probability  is  a 
strong  function  of  the  well  depth.  In  the  climate  models,  the  modulation  is  assumed 
to  result  from  a  weak  but  periodic  variation  in  the  eccentricity  of  the  earth’s  orbit 
with  period  ss  10’  yrs.  Digital  simulations  of  the  models  show  finite  amplitude  peaks 
in  the  power  spectra  located  at  the  modulation  frequency'*’. 

Several  years  ago,  a  physical  realization  of  stochastic  resonance  was  demon¬ 
strated  by  Fauve  and  Heslot®  using  an  electronic  Schmidtt  trigger  as  a  bistable 
system.  They  modulated  periodically  the  depths  of  the  potential  wells,  which  repre¬ 
sent  the  two  stable  states  of  the  switch,  alternately  and  simultaneously  applied  addi¬ 
tive  white  noise.  The  measured  power  spectrum  of  this  system  displayed  a  sharp 
peak  at  the  modulation  frequency  superimposed  on  a  slowly  varying  continuum  noise 
background  of  Lorentzian  line  shape.  These  authors  defined  the  signal-to-noise  ratio 
(SNR)  as  the  (measured)  amplitude  of  the  peak  relative  to  the  (measured)  noise  back- 
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ground.  They  observed  that  the  SNR  passed  through  a  maximum  as  the  noise  inten¬ 
sity  D  was  increased  from  zero,  and  that  at  the  maximum  the  value  D  -  Dg  could  be 
associated  with  the  period  T  of  the  modulating  function  by  a  Kramer’s  characteristic 
time  (see  also  Ref.  2),  t  »  T  oe  exp(Au/Do).  where  AU  is  the  unmodulated  barrier 
height. 

Very  recently,  McNamara,  Wiesenfeld  and  Roy  have  observed  this  phenomenon 
using  a  bidirectional  ring  laser  as  the  bistable  system’.  Using  an  intracavity  acousto¬ 
optic  modulator,  they  were  able  to  induce  changes  in  the  direction  of  the  lasing  by 
controlling  the  acoustic  frequency.  Representing  the  two  directions  as  the  stable  states 
separated  by  a  barrier,  it  was  possible  to  modulate  the  height  of  the  barrier  and  to 
introduce  external  noise  as  well  by  modulating  the  acoustic  frequency.  These  authors 
have  observed  the  same  phenomenology  as  those  of  Ref.  6.  In  particular,  they  have 
observed  a  sharp  peak  in  the  power  spectrum  of  the  laser  intensity  (measured  in  one 
dir''ction)  superimposed  on  a  broadband  noise  background  spectrum.  In  order  to 
obtain  the  SNR.  they  measured  the  amplitude  of  the  peak  and  that  of  the  noise  back¬ 
ground  at  the  modulation  frequency.  Measurements  of  the  SNR  versus  the  noise  in¬ 
tensity  in  this  experiment  demonstrated  the  characteristic  maximum,  though  the 
Kramer’s  time  was  not  obtained. 

In  this  Comment,  we  point  out.  based  on  physical  arguments,  that  the  power 
spectra  of  all  such  systems  regardless  of  the  details  of  the  model,  but  additively  mod¬ 
ulated  by  a  single  frequency  must  contain  delta  functions.  Noise  in  the  system  does 
not  alter  this  so  long  as  the  noise  does  not  multiply  the  amplitude  of  the  periodic 
modulation.  By  contrast,  experimental  measurements  of  the  power  spectra  or  digital 
simulations  of  models,  both  obtained  from  Fourier  transforms  of  time  series  made  at 
finite  resolution,  show  non-infinite  amplitude  peaks.  TTie  amplitudes  so  obtained  ref- 


-  3  - 


lect  both  the  strength  of  the  singularity,  which  contains  all  the  information  about  sto¬ 
chastic  resonance,  and  the  finite  resolution  of  the  Fourier  transform  technique  used. 
The  latter  is  an  entirely  instrumental  effect  which  also  determines  the  observed  line 
shape.  Results  obtained  by  measurements  of  the  amplitudes  only  of  the  peaks  in  the 
power  spectra  can  therefore  only  be  qualitative  in  the  absence  of  the  details  of  the 
line  shape®.  This  appears  to  have  escaped  notice  not  only  in  the  recent  experiments®-^ 
but  also  in  the  previous  digital  simulations®. 

For  the  purposes  of  discussion,  we  consider  two  models  represented  by  the  fol¬ 
lowing  Langevin  equations; 

X  -  ax  -  X®  +  fcoscjt  +  $(t).  £  <  (4a®/27)*/®.  (1) 


and 

X  -  a(t)x  -  X®  +  |(t).  a(t)  -  2(1  +  ecoswt)*/*,  e  <  1.  (2) 

which  represent  an  infinitely  damped  system  moving  in  the  standard  quartic  potential, 

U(x)  -  -(a/2)x®  +  (l/4)x®.  (3) 

modulated  at  the  frequency  u  and  driven  by  the  additive  noise  {(t).  Equation  (1)  rep¬ 
resents  additive  modulation  and  additive  noise.  In  this  case,  the  modulation  raises  and 
lowers  the  depth  of  each  well  alternately  on  alternate  half  cycles.  Equation  (2)  rep¬ 
resents  the  case  of  multiplicative  modulation  and  additive  noise®.  In  this  case,  the 
height  of  the  potential  barrier  AU  -  (l/4)a®  is  modulated  about  an  average  value  of 
1.  The  limits  on  e  insure  that  the  barrier  never  vanishes.  (In  the  simulation  discussed 
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below,  most  of  the  data  were  taken  for  a  -  (2)*/*  and  e  -  0.4). 

As  Eckmann  and  Thomas  have  pointed  out.  calculations  of  the  statistical  proper¬ 
ties  of  such  time  modulated  bistable  systems  are  by  no  means  trivial*®,  in  the  first 
instance  because  they  are  nonlinear,  and  in  the  second  because  they  are  not  station¬ 
ary.  Nevertheless,  some  general  remarks  can  be  made.  Considering  Eq.  (1).  whatever 
the  Fourier  transform  of  x(t)  may  be,  we  note  that  the  modulation  term,  ccoswt.  stands 
alone  and  therefore  will  necessarily  contribute  a  delta  function  at  u  to  the  transform 
and  hence  to  the  power  spectrum.  In  Eq.  (2),  or  in  other  versions  where  the  modula¬ 
tion  may  multiply  more  highly  nonlinear  terms,  the  situation  becomes  more  compli¬ 
cated  because  x(t)  itself  is  a  stochastic  function.  The  Fourier  transform  can.  in  princi¬ 
ple,  be  broadened  into  a  continuum  which  may  be  repeated  at  harmonics  of  the  mod¬ 
ulation  frequency.  In  the  simulation  of  Eq.  (2)  described  below,  however,  the  meas¬ 
ured  line  shapes  are  as  sharp  as  or  sharper  than  those  observed  in  the  simulation  of 
Eq.  (1).  In  ^either  case  is  the  line  shape  observed  to  be  continuously  broadened 
either  by  the  noise  or  by  the  dynamics  (see  Ref.  9).  Theoretical  results  recently  obta¬ 
ined  by  Fox  predict  Lorentzian  line  shapes  even  for  Eq.  (1)".  By  contrast,  a  general¬ 
ized  two  state  model  put  forth  by  McNamara  and  Wiesenfeld  predicts  delta  functions 
in  the  power  spectrum.** 

TTie  climate  models*"*  are  based  on  a  Langevin  equation  of  the  type  of  Eq.  (2). 
The  Schmidtt  trigger*  obeys  an  equation  of  the  type  of  Eq.  (1).  An  accurate  model 
for  the  laser  is  more  difficult  for  two  reasons.  First,  the  laser  is  a  multidimensional 
system,  and  even  with  appropriate  adiabatic  eliminations  probably  cannot  be  ade¬ 
quately  represented  by  any  less  than  two  coupled  Langevin  equations.  Second,  the 
details  of  the  bistable  potential  actually  introduced  by  the  AOM  are  not  well  under¬ 
stood**.  Nevertheless,  the  modulation  and  the  noise  must  enter  the  experimental  sys- 
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terns  phenomenologically  either  as  additive  or  multiplicative  terms  or  some  combina¬ 
tion  of  the  two.  We  believe  that  the  generic  systems  represented  by  Eqs.  (1)  and  (2) 
above  will  thus  reproduce  the  phenomenology  observed  in  the  experimental  systems. 

In  order  to  illustrate  these  remarks,  we  have  built  analog  simulators  of  Eqs.  (1) 
and  (2)  following  well  developed  techniques'*.  The  modulation  frequency  was  always 
set  at  f  -  cj/2ff  -  500  Hz.  and  the  modulation  amplitude  was  always  e  -  0.4.  An 
example  measured  time  series  x(t).  obtained  from  the  simulator  of  Eq.  (1)  is  shown  in 
Fig.  1  (a).  The  time  scale  is  delineated  by  the  plot  of  the  500  Hz  modulation  shown 
in  Fig.  1  (b).  The  switching  events  occur  randomly  in  time  and,  in  this  example,  on  a 
somewhat  longer  time  scale  than  the  modulation  period.  They  are.  however,  correlated 
with  the  modulation  as  shown  by  the  sharp  peak  at  exactly  500  Hz  in  the  measured 
power  spectrum  shown  in  Fig.  I  (c).  Figure  1  (d)  shows  a  power  spectrum  obtained 
from  a  simulator  of  Eq.  (2).  It  is  very  similar  except  that  a  small  peak  appears  at  the 
second  harmonic  of  the  modulation  frequency.  These  power  spectra  are  qualitatively 
very  similar  to  the  ones  published  by  Fauve  and  Heslol®  and  by  McNamara,  et.  al? 
in  the  sense  that  very  sharp  peaks  are  observed  superimposed  on  a  broadband  noise 
background. 

The  power  spectra  were  measured  in  the  following  way;  First  2048  points  of 
time  domain  data  x(t)  were  digitized  with  12  bit  accuracy,  each  point  separated  from 
its  neighbors  by  300  pis.  The  magnitude  of  the  square  of  the  Fourier  transform  was 
then  computed  and  compressed  into  1024  points.  The  final  power  spectrum  was  obta¬ 
ined  by  averaging  a  number  of  samples  (usually  200)  of  the  individual  spectra.  In 
every  case,  the  peak  at  the  modulation  frequency  was  only  a  few  points  wide,  and 
the  measured  amplitude  of  the  peak  depended  upon  which  individual  point  was 
chosen  as  "the  maximum".  Increasing  the  frequency  resolution  by  increasing  the  sep- 
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aration  time  between  digitized  points  resulted  in  narrower,  higher  amplitude  peaks 
and  the  reverse  was  also  true  though  we  explored  a  range  of  only  a  factor  of  two  in 
separation  time. 

Using  the  definition  adopted  in  Refs.  6  and  7  (the  ratio  of  the  amplitudes  of  the 
peaks  to  the  background  noise  power  density)  we  have  measured  the  SNR  as  a  func¬ 
tion  of  noise  intensity  D  -  where  t  is  the  noise  correlation  time'*,  for  the  sim¬ 

ulator  of  Eq.  (1).  The  data  are  shown  by  the  experimental  points  in  Fig.  2  and  are 
self  consistent,  since  all  points  were  measured  with  the  same  instrumental  resolution. 
However,  it  is  important  to  realize  that  lacking  a  quantitatively  measure  of  the  effect 
of  this  resolution,  i.e.  tacking  detailed  knowledge  of  the  line  shape,  the  vertical  scale 
in  Fig.  2  has  no  quantitative  meaning.  The  error  bars  shown  in  Fig.  2  require  dis¬ 
cussion.  They  were  assigned  by  making  repeatability  measurements  as  well  as  by 
testing  the  results  of  repositioning  the  cursors  presumably  located  at  the  maximum  of 
the  peak  and  near  the  base  of  the  peak  at  a  place  which  we  hoped  would  represent 
the  "average"  noise  power  density.  The  largest  amount  of  scatter  by  far  was  incurred 
with  the  repeatability  measurements.  The  reason  is  that  the  modulation  peaks  are 
only  a  few  points  wide,  so  that  small  variabilities  from  sample-to-sample  in.  for 
example  the  signal  generator  (modulation)  frequency,  resulted  in  large  variations  in 
the  peak  amplitudes  as  the  modulation  power  density  happened  to  be  shared  among  a 
few  or  several  "bins"  as  the  Fourier  transform  was  processed.  We  suggest  that  the 
comparably  large  error  bars  in  the  ring  laser  experimenf  might  have  derived  from 
the  same  variability. 

In  their  original  paper.  McNamara,  et  aV  outlined  a  theory  which  has  now  been 
written  up  in  detail'*.  We  quote  here  only  the  result; 
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S/N  -  (c/D2)exp(-2AU/D).  (4) 

where  S/N  is  the  power  density  amplitude  ratio  and  c  is  some  constant.  Following 
the  usual  definition,  the  SNR  in  decibels  (db)  is  given  by  SNR  -  101og(S/N).  As  did 
the  authors  of  Ref.  7.  we  have  found  it  necessary  to  add  an  offset  value  Dq  -  0.032 
vVHz  to  all  values  of  D  applied  to  the  simulator  in  order  to  fit  the  data  with  Eq.  (4). 
Presumably,  this  represents  the  effect  of  the  internal  circuit  noise.  For  c  -  60.8  and 
AU  -  0.5  the  theory  is  represented  by  the  curve  shown  in  Fig.  2  which  is  qualita¬ 
tively  comparable  to  that  obtained  in  the  Schmidlt  trigger  and  ring  laser  experiments. 

We  conclude  by  emphasizing  that  ideally  the  power  spectra  of  stochastic  reso¬ 
nance  systems  with  additive  modulation  contain  delta  functions  which  are  rendered 
into  measurable  peaks  by  experimental  systems  with  finite  frequency  resolution. 
Quantitative  measurements  of  the  SNR  of  such  systems  can  be  obtained  by  introduc¬ 
ing  the  details  of  the  instrumentally  broadened  line  shape,  or  alternatively  by  integ¬ 
rating  the  measured  power  spectra  thus  transforming  the  peaks  into  steps  whose  amp¬ 
litudes  are  independent  of  the  line  shapes. 

We  are  grateful  to  Kurt  Wiesenfeld  for  a  valuable  discussion.  Thanks  are  also 
due  to  Peter  Hanggi  for  remarks  concerning  the  Four’er  transforms  of  modulated  sys¬ 
tems.  This  work  was  supported  by  the  Office  of  Naval  Research  grant  no. 
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FIGURE  CAPTIONS 


Fig.  1.  Example  results  of  the  analog  simulation,  (a)  A  time  series  for  Eq.  (1)  with  a 
-  (2)^/*,  D  -  0.20  v*/H2.  and  t  -  0.20.  (b)  The  500  Hz  modulation  which  also  estab¬ 
lishes  the  time  scale  for  (a)  with  €  -  0.4.  (c)  An  example  power  spectrum  which  is 
the  result  of  200  averages  computed  from  individual  time  series  obtained  for  the  con¬ 
ditions  stated  in  (a)  and  (b).  TTie  sharp  peak  at  500  Hz  establishes  the  frequency 
scale,  (d)  A  power  spectrum  obtained  for  the  same  conditions  as  listed  in  (a)  except 
from  the  simulator  of  Eq.  (2),  and  only  100  averages  were  accumulated. 

Fig.  2.  The  measured  SNR  versus  D‘A  from  the  simulator  of  Eq.  (I)  with  a  -  (2)*/* 
(which  implies  that  AU  -  0.5),  t  «■  0.20  and  the  modulation  as  in  Fig.  1  (b).  For  a 
discussion  of  the  error  bars,  see  the  text.  The  curve  is  Eq.  (4)  with  c  -  60.8  and  AU 
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ABSTRACT 

We  have  constructed  an  electronic  simulator  of  a  simple  bistable  system 
driven  by  noise,  whose  intensity  is  determined  by  the  instantaneous  value  of 
the  coordinate.  We  observe  that  the  most  probable  state  of  the  system  can  be 
reversed  by  altering  the  noise  intensity  only  in  the  neighborhood  of  the 
barrier;  an  effect  pointed  out  by  Landaiuer  many  years  ago  in  the  context  of 
discussions  on  entropy  related  stability  criteria  for  nonequilibrium  systems. 
We  compare  detailed  measurements  on  the  system  with  the  recent  white  noise 
calculations  of  Landauer  and  van  Kampen.  The  system  also  has  interesting 
possibilities  for  tests  of  contemporary  colored  noise  theory  which  we 
illustrate  with  an  example. 

PACS  05.40.+j 


1 


I .  Introduction 

Many  years  ago  Landauer  raised  a  discussion  on  the  implications  of  state 
dependent  diffusion  in  the  context  of  entropy  generation  and  stability  in 
nonequilibrium  systemsJ  His  point  was  that  the  relative  stability  of  a 
multistable  system  could  be  altered  by  path  dependent  diffusion  even  when 
localized  to  the  neighborhoods  of  the  potential  barriers  separating  the 
deterministically  stable  or  metastable  modes.  That  is,  the  most  probable 
state  of  such  a  system  cannot  be  determined  from  information  about  the  local 
minima  alone.  Later  it  was  observed  that  in  certain  models,  state  dependent 
diffusion  can  even  generate  maxima  in  the  probability  distribution  at 
locations  where  no  potential  minima  exist.2  These  early  studies  gave  rise  to 
wide  discussions  which  continue  to  be  of  current  interest. ^ 

The  more  general  problem  of  diffusion  in  inhomogeneous  media  has  recently 
been  carefully  analyzed  by  van  Kampen.^  State  dependent  diffusion  which  is 
.periodic  and  commensurate  with  a  spatially  periodic  potential  can  give  rise  to 
continuous  currents  in  the  absence  of  externally  applied  fields.  Such  systems 
have  recently  been  treated  by  Buttiker^  and  van  Kampen^  as  well  as  Landauer.^ 
In  this  paper  we  consider  the  simplest  model  for  state  dependent 
diffision  in  a  bistable  system  as  originally  suggested.^  We  take  the  standard 
quart! c  as  the  potential 

V(x)  *  -  x2/2  +  x^/4  +  EX,  (la) 

with  the  dynamics  being  determined  by  the  infinitely  damped  system 

X  “  X  -  x^  +  g(x)c(t)  -  E,  (lb) 

where  C(t)  is  a  Gaussian,  white  noise  with  zero  mean  and  correlation 
<C(t)5(s)>  *  2nn6(t-s);  and  where  g(x)  is  defined  by 
g(x)  »  1,  for  X  <  0  or  X  >  1/2 

g(x)  =  1  +  A,  for  0  <  X  <  1/2.  (2) 
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The  noise  intensity  is  and  e  controls  the  symmetry  of  the  potential  wells. 
This  specific  system  has  recently  been  analyzed  by  both  van  Kampen,®  who,  in  a 
notably  clear  tutorial,  derived  the  probability  densities  from  thermodynamic 
arguments  as  well  as  from  the  diffusion  equation,  and  by  Landauer.^ 

The  effect  of  the  reversal  of  stability  in  this  bistable  system  due  to  an 
increased  temperature  (larger  D)  in  the  neighborhood  of  the  barrier  is 
illustrated  in  Fig.  1,  Though  the  illustration  is  schematic,  the  graphs  are 
quantitative.  On  the  left  is  shown  U(x)  as  given  by  Eq.  (la)  with  e  =  0.05. 

On  the  right  are  shown  the  probability  densities  as  measured  on  the  simulator 
described  in  Section  III  below  for  Increasing  values  of  a.  The  effect  of  A  in 
reversing  the  most  probable  state  is  clearly  evident. 

This  paper  is  organized  as  follows:  In  Section  II  the  theory  is  briefly 
reviewed.  Expressions  for  the  ratio  of  the  amplitudes  of  the  stationary 
probability  densities  in  the  two  wells  and  for  the  magnitudes  of  the  discon¬ 
tinuities  at  the  locations  where  the  noise  intensity  Is  discontinuously  changed 
are  summarized.  In  Section  III  the  simulator  and  the  methods  of  measurement 
are  described.  The  main  results  are  presented  in  Section  IV  where  measurements 
of  the  amplitude  ratio  of  the  probability  density  and  the  magnitudes  of  the 
discontinuities  are  compared  to  the  theory  of  Refs.  6  and  7,  In  Section  V  we 
examine  the  colored  noise  problem  and  the  possibility  that  this  system  could  be 
of  use  in  the  testing  of  contemporary  colored  noise  theory.  As  examples,  data 
on  the  amplitude  ratios  are  compared  to  the  predictions  of  the  conventional 
small  correlation  time  theory,  as  recently  improved  by  Fox,^  and  to  Hanggi's 
ansatz The  advantage  of  this  system  is  that  the  results  depend  only  on  the 
ratios  of  probability  densities,  so  that  the  correlation  time  dependent 
prefactors  which  often  appear  in  various  forms  in  colored  noise  approximate 
theories  cancel  out,  revealing  the  exponential  behavior  alone. 
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This  could  be  a  considerable  advantage,  since  it  is  often  difficult  or 
impossible  to  clearly  distinguish  between  the  influence  of  the  prefactors  and 
exponentials  using  numerical,  matrix  continued  fraction  or  analog  simulations. 
A  vigorous  discussion  on  the  merits  and  accuracy  of  various  colored  noise 
approximate  theories  is  currently  in  progress. Finally  in  Section  VI  we 
summarize  our  results. 

II .  Theory 

We  consider  the  symmetric  potential  defined  by  Eq.  (la)  with  e  =  0  and  as 
shown  in  Fig  (2b).  Following  Landauer's  notation, ^  the  left  hand  well  is 
located  at  A  (x  =  -1),  the  barrier  at  8  (x  =  0),  the  region  of  increased  noise 
intensity  between  B  and  C  (x  =  1/2)  and  the  right  hand  well  at  D  (x  =  1). 

It  is  necessary  to  first  consider  what  happens  at  a  temperature  (or  noise 
intensity)  discontinuity.  For  stationary  conditions  the  probability  currents 
across  the  discontinuity  at  B,  for  example,  can  be  written  p(B")v(B")  * 
p(B“*’)v(B'^) ,  where  p  is  the  probability  density  and  v  is  the  velocity  of  a 
particle  in  thermal  equilibrium  at  temperature  T  just  to  the  left  of  the 
discontinuity  (B-)  or  just  to  the  right  (B'*’).  The  velocity  is  v  ,  /~T, 
so  that  p(B‘*’)/p(B")  =  (Tl/Th)^/2,  where  Tl  is  the  (lower)  temperature  outside 
the  region  BC  and  Th  is  the  (higher)  temperature  inside  BC.  We  identify  the 
temperature  with  the  diffusion  or  noise  intensity  through  n^,  »  pkT,  and  since 
we  consider  only  homogeneous  media  we  take  p  =  const.  =  1.  The  ratios  are  thus 
p(B+)/p(B-)  »  p(C-)/p(C+)  =  (OnL/nnH)^/^  (3) 

Landauer  then  argues  that  the  densities  within  the  regions  are  given  by 
th’  Boltzmann  distribution  exp(-U(x)/nf,) .  In  addition  to  Eq.  (3),  the  ratios 
are  p(B")/p(A)  =  expC-(UB-U;y)/npL];  p(C.)/p(B''‘)  =  expC-(Uc-IJB)/nnn];  ^nd 
p(D)/p(C‘*')  a  expC-(Uo-U(.)/OnL].  These  probabilities  are  then  multiplied 
together,  whereupon  the  ratios  given  by  Eq,  (3)  cancel,  and  the  result  is 
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p(D)/p(A)  =  expC-(Uo-UA)/DnL]  expC-AU(DnL-DnH)/OnLnf,^] , 
where  aU  =  Uq  -  Ur.  van  Kanpen^  obtains  the  same  result  for  the  amplitude 
ratio  at  the  wells,  but  he  predicts  that  for  ratio  at  the  discontinuities 
p(B+)/p(B-)  »  p(C-)/p(C+)  -  DnL/OnH  (5) 

instead  of  Eq.  (3).  In  this  work,  we  consider  only  the  symmetric  (e  *  n) 
potential.  With  =  Uq,  Eq.  (4)  becomes 

p(D)/p(A)  =  exp[-AU(DnL-DnH)/*^nL®nH^* 

III.  The  Simulator 

Figure  2(a)  shows  a  schematic  diagram  of  the  simulator  of  Eqs.  (1)  and 
(2).  The  design  is  straight-forward  and  one  which  we  have  used  before 
with  the  exception  of  the  added  system  for  generating  g(x).  This  is  shown  by 
the  two  comparitors  which  continuously  test  the  voltage  on  x  and  compare  it  to 
the  preset  values  Xq  and  x-;  which  mark  the  boundries.of  the  region  of 
increased  noise  intensity.  Each  time  the  trajectory  crosses  one  of  the 
boundries  an  output  is  applied  to  the  AND  gate  shown.  This  gate,  in  turn,  and 
the  adder  which  follows  it,  provide  a  suitable  voltage  Vg,  to  a  multiplier 
whose  other  input  is  the  noise  voltage  V^.  The  logic  is  arranged  in  such  a 
way  that  Vg  =•  1  +  A  when  Xq  <  x(t)  <  X] ,  and  Vg  =  1  when  x(t)  is  outside  this 
range,  or  Vg  =  g(x).  The  noise  voltage  lit)  which  drives  the  trajectory  in 
the  bistable  potential  is  then  defined  by 

5(t,x)  »  g(x)Vn(t)  ,  (7) 

with  g(x)  defined  by  Eqs.  (2). 

The  noise  voltage  is  supplied  by  a  noise  generator^^  of  wide  hut  finite 
bandwidth,  and  is  therefore  necessarily  colored.  In  order  to  define  its 
correlation  time  tq  precisely,  it  is  passed  through  a  linear,  single  pole 
filter  with  transfer  function  H(a))  »«  1/Cl  +  The  simulator  scales 
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time  with  the  integrator  time  constant  t-f ,  so  that  the  dimensionless 
correlation  time  is  t  =  Tn/f-j ,  and  the  correlation  function  of  is 

<Vn{t)Vn(s)>  =  (Dn/T)exp(-|t  -  s|/t)  (ft) 

In  the  limit  t  -*■  s,  the  noise  intensity  is  defined  by 

On  =  <  V  >.  *  (9) 

and  the  mean  square  noise  voltage  is  a  measured  quantity.  In  this  simulation 
=  ino  ^is  and  the  range  of  was  20  ps  to  500  jxs  (the  range  of  x  was  then 
0.20  to  5).  For  comparison  with  the  white  noise  theories,  however,  we  set  the 
value  X  =  0.20.  Equations  (2),  (7)  and  (9)  then  give  for  the  noise  intensities 


f^nL  ’  ^  <  ''^2  > 

OnH  =  -c  <  Vn2  >  (i  +  a)2^  (10) 

Figure  2(b)  shows  the  potential  of  Eq.  (la)  with  e  =  0.  The  region  RC  of 
enhanced  noise  intensity  is  shown  by  the  darkened  segment  which  lies  between 
Xq  =  0  and  xi  =  1/2.  The  potentials  shown  in  Eq.  (4)  are  then  given  by  * 

-  1/4,  Ub  =  0,  Dq  »  7/64  and  Uq  =  1/4,  as  determined  hy  Eq.  (la)  with  e  =  0. 
The  theoretical  prediction  for  the  amplitude  ratio  of  the  densities  in  the 
wells,  Eq.  (4),  then  becomes 


7 

p(0)/p(A)  =  exp  I - - 


(11) 


Using  Eqs.  (3),  and  (10)  the  predictions  for  the  amplitude  ratios  at  the 
discontinuities  are 


1/(1+A) 

p(B+)/p(B-)  »  p(C-)/p(C'')  =  or 

1/(1+A)2  (12b) 


The  simulator  is  operated  by  first  setting  a  =  0  and  hy  adjusting  e  in 
order  to  obtain  a  symmetric  density  p.  Checks  on  this  symmetry  adjustment 
were  frequently  made  throughout  the  experiment.  Figure  3  shows  an  example 
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set  of  measured  densities  for  increasing  values  of  a  commencing  with  the 
symmetric  density  shown  at  the  top.  These  densities  were  assembled  as 
averages  from  a  sequence  of  time  series  x{t)  of  4096  digitized  points  each. 
Typically  ROO  such  time  series  were  obtained  so  that  each  density  resulted 
from  about  3.2  x  10®  digitized  points.  For  each  density,  the  amplitudes  at 
the  two  peaks  and  at  the  high  and  low  sides  of  each  discontinuity  were 
recorded. 

IV.  The  results 

Figure  4  shows  the  results  of  our  measurements  of  the  amplitude  ratio  in 
the  wells  as  a  function  of  A  for  various  values  of  D^.  The  solid  lines  are 
plots  of  Eq.  (11)  with  OpL  determined  directly  from  x  and  measurements  of  < 

>.  There  are  no  adjustable  constants.  The  agreement  between  Eq.  (11)  and 
our  measurements  is  excellent,  but  probably  fortuitous,  since  the  noise 
driving  our  simulator  Is  actually  colored  with  t  ®  0.2,  and  the  theory  is 
valid  only  for  white  noise.  The  region  of  quasi-white  behavior  of  the 
simulator  is  x  «  1,  but  in  practice  the  usable  dynamic  range  of  the  analog 
components  imposes  a  lower  limit  of  x  -  0,  1  to  0.2  for  reliable  operation 
without  voltage  clipping.  Nevertheless  the  systematic  behavior  with  D^l  and  a 
is  convincing. 

The  scatter  in  the  data  points  shown  on  Fig.  4  is  the  result  of  both 
systematic  and  statistical  errors.  The  ratio  p(0)/p(A)  is  obtained  from  two 
large  numbers,  each  with  a  certain  statistical  error.  The  statistical  scatter 
was  estimated  from  short  term  repeatability  measurements,  and  is  shown  as  the 
example  error  bar  on  the  =  0.36  data.  A  more  troublesome  error  stems  from 
longer  term  variations  in  <  >  due  to  very  low  frequency  drifts,  or 

systematic  variations,  in  the  noise  generator.  An  attempt  has  been  made  to 
estimate  the  maximal  effects  of  this  error,  and  is  shown  by  the  dashed  curves 
for  OpL  *  0,40  (upper  curve)  and  D^l  “  0.32  (lower  curve). 
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We  have  also  measured  the  amplitudes  of  the  density  on  both  sides  of  the 
discontinuities  in  order  to  find  the  ratios  p(B''')/p(B")  and  p(C")/p{C'’') .  Both 
Refs.  (6)  and  (7)  make  specific  predictions  for  the  behavior  of  these  ratios 
with  A  as  shown  by  Eqs.  (12).  In  Ref.  (7)  it  is  argued  additionally  that 
whatever  fhe  functional  dependence  of  the  magnitude  of  the  discontinuities  on 
A,  it  can  be  expected  to  be  the  same  for  both  discontinuities,  and  therefore 
will  cancel  in  the  result  Eq.  (11).  Our  measurements  of  these  ratios  are 
shown  on  Fig.  5,  where  (a)  and  (b)  display  the  results  for  a  larger  and  a 
smaller  noise  intensity  respectively.  The  statistical  errors  on  these  data 
are  relatively  larger,  as  can  he  judged  from  the  scatter,  because  the 
discontinuities  occur  at  low  amplitude  on  the  density.  The  theoretical 
predictions  are  shown  by  the  curves;  Ref.  (7)  and  Eq.  (12a)  are  the  solid 
curves,  and  Ref.  (6)  and  Eq.  (12b)  are  the  dashed  curves. 

We  remark  that  both  theories  predict  the  equality  of  the  magnitudes  of 
the  two  discontinuities,  and  their  independence  on  0.  Certainly  to  within  the 
scatter  of  our  data  the  measured  magnitudes  appear  to  be  equal  as  shown  by  the 
solid  circles  (discontinuity  at  R,  x  =  0)  and  the  open  circles  (discontinuity 
at  C,  X  =  1/2).  Both  sets  of  data  seem  to  lie  systematically  higher  than  Eq. 
(12a),  a  trend  which  seems  more  evident  for  large  noise  intensity. 
Nevertheless,  the  systematics  clearly  favor  Eq,  (12a), 

Certainly  these  results  come  as  no  surprise  to  anyone,  since  exact 
solutions  of  white  noise  systems  in  the  limit  of  large  damping  are  very  well 
known.  The  quasi-white  noise  measurements  shown  here  serve  to  indicate  the 
accuracy  of  our  simulator. 
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V.  Colored  Noise 

Nearly  all  colored  noise  approximate  theories  begin  with  an  effective 
Fokker-PI anck  equation 

2 

—  p(x,t)  =  [IJ' (x)p(x,t)]  +  D(x)p(x,t),  (13) 

at  ax  ax2 

wherein  the  diffusion  n(x)  nfx.T),  is  intended  to  approximately  account  for 

the  nonzero  point-to-point  correlations  over  x  which  are  induced  by  a 

correlated  driving  force  (the  colored  noise).  [The  case  wherein  n(x)  is  the 

result  of  an  inhomogeneous  medium  has  been  recently  examined  by  van  Kampen^l. 

The  recent  rapid  growth  of  the  theoretical  colored  noise  industry^"^^  has 

resulted  in  a  variety  of  expressions  for  the  "renormalized"  diffusion  n(x,-c). 

For  the  purpose  of  example  we  choose  here  only  two:  the  improved  "small 

t"  approximation  of  Fox®>^^  which  reduces  to  the  often  cited  result  of  Sancho, 


et  £l_^5  in  the  limit  of  small  t;  and  the  ansatz  due  to  Hanggi  The  Fox 
results  are  obtained  from  Eq.  (13)  with 

Df(x,'c)  =  DgCl  +  T  lJ"(x)]-^  ,  ■  (14) 

which  results  in  the  stationary  density 

pp(x)  =  (1  -  T  +  3x2x)exp[-  UF(x,T)/nQ]  ,  (15a) 

with  Uf(x,'c)  for  our  potential  given  by 

Uf(x.t)  =  -  x2/2  +  x4/4  +  -c  (x2/2  -  x^  +  x®/2)  (15h) 

Hanggi's  results  are  obtained  from  Eq.  (13)  with 

Dh(x,t)  =  noCl  +  'c(3<  x2  >  -  1)]-1  ,  (Ifi) 

where  <  x2  >  is  an  average  over  the  trajectories,  and  this  n(x)  is  specific  to 
our  potential.  The  density  in  this  case  is 

PH(x)  =■  N  expC-l)(x)/DH(x,'c)]  (17) 

where  U(x)  is  given  by  Eq.  (la). 
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As  a  zeroth  order  approximation,  we  can  simply  substitute  these  results 
into  the  white  noise  formula  for  the  ratio  of  the  densities.  This  procedure 
is  based  on  the  observation  that,  even  though  this  is  a  system  with  state 
dependent  noise,  the  space  is  separated  into  three  regions,  within  each  of 
which  the  noise  intensity  is  constant.  Nevertheless,  to  simply  substitute 
colored  noise  densities  for  the  white  ones  neglects  the  effects  (if  any)  of 
correlations  across  the  boundries.  Our  only  justification  for  this  procedure 
is  that  our  measurements  of  the  amplitude  ratios  of  the  discontinuities,  as 
discussed  below,  show  no  systematic  behavior  with  t  to  within  our  (not  so 
small)  statistical  errors. 

Using  Fqs.  (15)-(17),  results  in  the  following  predictions 


p(D)/p(A)p  *  exp((-- — )[1  - 

OnL 


1 

(777)2 


]} 


for  the  Fox  theory;  and 


(18) 


p(D)/p(A)h  =  exp{- 


7(1  2t) 

640nL 


[  1 


1 

(1  +  a)^ 


]} 


(19) 


for  the  Hanggi  theory.  We  note  that  the  t  dependent  prefactors  (if  any)  on 
the  densities  cancel  out  in  this  application.  Further,  these  two  expressions 
make  qualitatively  opposite  predictions:  Eq.  (18)  shows  the  ratio  increasing 
with  T,  while  Eq.  (19)  shows  it  decreasing.  Both  show  the  same  A  dependence 
(which  is  also  the  white  noise  A  dependence). 

Our  measurements  at  fixed  a  =  1  are  shown  in  Fig.  fi  where  the  open 
circles  are  for  the  large  O^l  and  the  closed  circles  are  for  the  small  np^. 
Equation  (18)  is  shown  by  the  dashed  curve  and  Fq.  (10)  by  the  solid  curve. 
Both  are  plotted  for  Op^  =  0.40  which  is  also  the  noise  intensity  for  which 
the  solid  circles  were  obtained.  The  lozenge  on  the  vertical  axis  shows  the 
white  noise  limit  for  all  theories  and  illustrates  the  difficulty  in 
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* 

comparing  our  simulations  with  white  noise  results:  the  quasi-white  noise 
simulation  was  done  for  t  =  (1.2  but  does  not  extrapolate  well  to  the  lozenge 
shown  on  the  axis  at  -c  =  0. 

Finally,  measurements  of  the  amplitude  ratios  at  the  discontinuities 
versus  x  for  A  =  1  are  shown  in  Fig.  7.  There  is  no  discernible  x  dependence 
of  these  ratios. 

VI .  Summary  and  conclusions 

We  have  measured  stationary  density  amplitude  ratios  for  a  system  with 
state  dependent  noise  applied  in  d  of  three  discrete  regions.  The  noise 
intensity  was  changed  discontinuously  at  the  boundries  of  the  inner  region. 

The  measured  results  are  in  good  agreement  with  the  predictions  of  white  noise 
calculations  due  to  Landauer  and  van  Kampen.  In  addition,  we  have  repeated 
two  sets  of  amplitude  ratio  measurements  for  a  range  of  noise  correlation 
times,  and  compared  the  results  to  two  current  colored  noise  approximate 
theories.  While  neither  approximation  accurately  describes  the  data,  the 
ansatz  of  Ref.  9  is  in  better  qualitative  agreement. 

We  are  grateful  to  R.  Landauer  and  M.  Ruttiker  for  stimulating 
discussions  and  for  suggesting  this  experiment.  We  are  also  indebted  to  R. 
Fox,  P,  Hanggi  and  J.  M,  Sancho  for  various  discussions.  This  work  was 
supported  by  the  Office  of  Naval  Research,  Grant  No,  N00014-88-K-00R4. 
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FIGURE  CAPTIONS 


Fig.  1  Potentials  from  Eq.  (1)  with  e  =  0.05  are  shown  on  the  left  and 
measured  densities  on  the  right  for  A  =  0,  0.45,  and  1.3 
top-to-bottom,  showing  the  reversal  of  the  most  probable  state. 

Fig,  2  (a)  A  schematic  diagram  of  the  simulator,  showing  multipliers  (x) 

summers  (±),  comparltors  (comp.)  and  an  and  circuit  (AND),  (b)  An 
example  of  the  potential  utilized  in  the  simulator  with  the  region  of 
increased  noise  intensity  BC  shown  in  hold  lines  between  Xq  =  n  and 

XI  «  1/2. 

Fig.  3  Probability  densities  measured  for  values  of  a  =  0,  0.5,  0,75,  and 
1.0  top-to-bottom.  The  region  of  Increased  noise  Intensity  RC  Is 
clearly  evident  for  A  >  0. 

Fig.  4  p(0)/p(A)  versus  A  for  values  of  0^(_  »  0.20  (squares),  0.3fi 

(triangles),  O.fiO  (circles)  and  1,6  (Inverted  triangles).  The  solid 
lines  are  plots  of  Eq.  (11).  The  dashed  lines  estimate  the 
systematic  error  In  DnL» 

Fig,  5  p(B'*’)/p(B*)  (solid  circles)  and  p(C")/p(C+)  (open  circles)  for  (a) 

OpL  *  1*60  (b)  OpL  ®  0.36,  The  solid  line  is  Eq.  (12a)  and  the 

dashed  line  Is  Eq.  (12b). 

Fig.  6,  The  effect  of  colored  noise  Is  shown  by  this  plot  of  p(D)/p(A)  versus 

X  for  A  »  1.  The  open  circles  are  for  O^l  *  1*6  and  the  solid 

circles  are  for  O^l  *  0,40.  The  dashed  curve  is  Eq.  (18)  and  the 
solid  curve  Is  Eq.  (19),  both  for  OpL  *  0,40, 

Fig.  7.  The  effect  of  noise  color  x  on  the  discontinuities  p(B'*')/p(R") 

(solid  circles)  and  p(C“)/p(C*)  (open  circles).  The  upper  data  set 
Is  for  OpL  ■  1 .60  and  the  lower  set  for  D^l  ®  0.40,  All  data  are  for 
A  *  1 . 
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ABSTRACT 


We  show  via  an  approximate,  analytical  technique  and  an  analog  simulation  of 
the  classical  Langevin  equation  of  a  correlated,  spontaneous  emission  laser  (CEL)  that 
noise  of  non-zero  correlation  time  leads  to  an  enhancement  in  the  characteristic  CEL 
noise  quenching. 


PACS  42.50.-p;  42.65.-k;  05.40.+j 
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The  recognition  that  quantum  noise  determines  the  ultimate  accuracy  of  active 
interferometers  in  applications  involving,  for  example,  gravitational  wave  detection^"*, 
or  ring-laser  gyroscopes*'®,  has  stimulated  the  investigation  of  laser  systems  with 
reduced  spontaneous  emission  noise  such  as  the  correlated  (spontaneous)  emission 
laser®”®  (CEL).  In  the  CEL.  the  relative  phase  angle  between  two  electromagnetic 
waves  can  be  freed  completely  from  the  effects  of  spontaneous  emission  noise  as  a 
result  of  the  fact  that  this  noise  appears  to  the  phase  as  a  purely  multiplicative 
stochastic  process*.  The  physics  of  noise  quenching  in  the  CEL  has  been  discussed 
in  parts  I  and  II  of  this  series  of  papers*®.  Moreover,  it  has  been  shown  that,  in  a 
laser  consisting  of  long-lived  atoms,  the  quantum,  i.e.  the  spontaneous  emission, 
noise  is  colored  instead  of  white.**  An  excellent  discussion  of  both  white  and  col¬ 
ored  (pump)  noise  effects  in  lasers  has  been  provided  recently*®.  In  this  brief  note, 
we  demonstrate  the  purely  classical  quenching  effect  of  colored,  multiplicative  noise 
in  a  system  with  a  periodic  potential  as  has  been  shown  to  apply  to  the  CEL.  Such 
effects  have  been  previously  studied  theoretically  for  white,  additive  noise  using  con¬ 
tinued  fractions*®  and  for  white,  multiplicative  noise  by  approximate  analytic  tech¬ 
niques*®. 

We  do  not  intend  to  present  here  a  theory  of  the  CEL  with  long-lived  atomic 
states,  nor  are  we  concerned  with  the  intricacies  of  colored  noise  theory  and  the  var¬ 
ious  approximation  schemes  associated  therewith**  as  used  to  calculate  approximate 
Fokker-Planck  equations  (except  as  in  Ref.  22).  Rather,  we  focus  on  the  influence  of 
noise  color  on  the  noise  quenching  effect  as  obtained  from  the  classical  Langevin 
equation  of  the  two-mode  CEL.  A  similar  equation  also  describes  the  phase  of  the 
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electromagnetic  field  in  a  two-photon  CEL‘®. 

We  pursue  this  study  following  two  strategies:  First,  we  perform  an  approxi¬ 
mate  analysis  of  the  Langevin  equation  by  linearization  within  a  physically  relevant 
range  of  parameters.  With  the  help  of  5ome  standard  techniques  of  noise  theory**,  we 
find  the  (approximate)  width  of  the  steady  state  probability  distribution,  Po(^)  for  the 
relative  phase  0,  to  be  reduced  by  a  factor  (1  +  yr)~^  compared  to  the  white  noise 
case.  Here,  t  is  the  noise  correlation  time  and  7  -  (b*  -  a*)*/*,  where  b  and  a  denote 
the  laser  gain  and  the  detuning  between  the  two  waves  respectively.  Thus  the  CEL 
noise  quenching  is  not  only  preserved  in  the  case  of  colored  noise  but  is  even 
enhanced.  Moreover,  the  white  noise  results*®  for  the  noise  induced  drift  <A>,  and 
the  noise  induced  asymmetry  are  both  reduced  by  this  factor.  In  the  second 
approach,  we  simulate  the  Langevin  equation  of  the  CEL  with  an  electronic  circuit 
using  a  circuit  and  techniques  similar  to  those  previously  described*®..  The  simulator 
provides  immediately  the  steady  state  distributions  Pq  as  shown  by  the  examples  in 
Fig.  1,  for  a  range  of  values  of  t.  We  compare  and  contrast  these  distributions  to 
Gaussian  approximations  with  first  and  second  moments  <A*’*>.  We  find  qualitative 
agreement  between  the  measured  and  calculated  results  and,  in  particular,  we  are 
able  to  confirm  the  predicted  enhancement  of  the  CEL  noise  quenching  by  color. 

The  natural  extension  of  the  CEL  Langevin  equation  for  the  relative  phase 
difference  between  the  two  modes  in  the  presence  of  Gaussian  noise, 

<e(t)e(s)>  -  (D/T)exp[-lt  -  s|/t].  (1) 

of  noise  intensity  D  and  correlation  time  r  with  zero  mean. 
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<e(t)>  -  0. 


(2) 


islO.17,18 


0  -  a  +  bsin0  +  e(t)sin0.  (3) 

For  small  noise  intensity  D/b  «  I,  and  for  small  detuning  |a|/b  «  1.  we  can  try  the 
ansatz^^'^’’ 


0(t)  ff  +  arcsin(a/b)  +  A(t),  (4) 


where  |A|  «  1.  Thus  Eq.  (3)  reduces  to 

A  -  -[7  +  (7/b)€(t)]A  -  (a/b)€(t). 


where  7  was  defined  above  in  terms  of  a  and  b,  and  with  the  obvious  solution 


A(t) 


Aoexp 


-I 

b 


r  ^ 

dim 

Jo 


f  ‘ 

-t 

dt'c(t')exp 

Jo 

-Hi  -  tvj 

df'€(t") 

t' 

(5) 


where  Aq  -  A(t-O).  From  Eq.  (5)  it  is  straight  forward  to  evaluate  the  moments^®  (AJ) 
for  j  -  1  and  2. 
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With  the  help  of  Eqs  (1)  and  (2)  we  can  find  the  center-of-gravity.  <0>,  for  the 
stationary  distribution,  at 

<0>  s  )r  +  arcsin(a/b)  +  <A> 

S!  ff  +  arcsin(a/b)  +  (I  +  7T)"‘(a/b)(D/b).  (6) 

and  the  approximate  width  of  Pq  is  governed  by 

5  <02>  -  «f»‘^  a  (1  +  7T)~*(a/b)^(D/7).  (7) 

The  Gaussian, 

Po^PP  (0)  -  (2ir)-‘/2a-‘exp[-(2(ry‘(0  -  <0»^].  (8) 

with  <0>  and  a  given  by  Eqs.  (6)  and  (7),  thus  represents  the  simplest  approximation 
to  the  exact  distribution  Pq,  which  we  do  not  calculate  here. 

The  noise  source  €(t)  is  multiplicative®  and  so  gives  rise  to  a  noise-induced  asym¬ 
metry  which  manifests  itself  as  a  separation  of  the  center-of-gravity  of  the  right 
hand  peak  of  Pq.  given  by  Eq.  (6),  and  the  location  of  its  maximum^*, 

•^max  -  If  +  arcsin(a/b)  -  (1  +  7Tr‘(a/b)(D/b).  (9) 

This  asymmetry  is  obviously  not  contained  in  the  Gaussian  approximation. 


-  7  - 


From  Eqs.  (6),  (7)  and  (9).  we  recognize  that  the  noise  color  modifies  the  white 
noise  (t  -  0)  moments  by  the  prefactor  (1  +  yT)~K  which  signifies  an  enhancement  of 
the  noise  quenching  characteristic  of  the  CEL.  We  emphasize  that  this  noise  color 
enhanced  quenching  is  an  entirely  classical  effect  which  may  have  applications  to 
other  systems.  Nor  does  the  quenching  stem  from  the  multiplicative  nature  of  the 
noise  in  this  particular  application.  For  example,  a  similar  result  can  easily  be  obta¬ 
ined  for  the  standard  phase  locked  ring  laser  Langevin  equation^®, 

0  =  a  +  bsin0  +  c(t).  (10) 

For  this  case,  following  an  analysis  similar  to  the  above  results  in  an  approximate 
width  for  the  stationary  distribution  in  the  phase  locked  region  given  by, 

"ring  laser  ^  (1  *  Cl) 

which  shows  that  the  width,  for  this  simpler  case  of  additive  noise,  is  governed  by 
the  same  color  prefactor. 

The  techniques  for  simulating  periodic  potentials  with  analog  circuits  have  been 
previously  discussed*®,  and  here  we  use  the  same  approach  except  for  the  following 
simple  alterations:  First,  the  hybrid  analog-digital-analog  system  for  producing  the 
spatially  periodic  force  was  replaced  with  an  analog  chip^®  which  accomplishes  the 
same  task;  and  second,  an  additional  multiplier  was  used  to  implement  the  multiplica¬ 
tive  noise.  The  voltages  representing  0(t)  were  digitized  into  time  series  of  approxi¬ 
mately  four  million  points  from  which  the  stationary  probability  densities  were 
assembled.  In  contrast  to  the  simulation  of  this  same  system  with  additive  noise 
described  in  Ref.  16,  in  the  present  case  only  the  one-dimensional  densities  were 


-  8  - 


measured.  Our  purpose  is  simply  to  demonstrate  the  phenomenon  of  noise  quenching 
and  the  noise  color  effects  in  this  purely  classical  multiplicative  system. 

Figure  1  shows  a  set  of  three  probability  densities  measured  for  values  of  t 
ranging  from  0.1  (quasi-white  noise)  to  10  (strongly  colored  noise).  These  densities  are 
consistentconsistant  with  the  behavior  predicted  by  Eqs.  (6)  and  (9).  which  predict 
that  the  effect  of  the  multiplicative  noise  is  to  induce  a  r-dependent  shift  between 
the  center-of-gravity  and  the  maxima  of  the  peaks  of  Pq.  This  shift  is,  however, 
quite  small  for  small  D  and  a  where  the  theory  should  be  accurate.  For  the  values  of 
the  parameters  shown  in  the  Figure  caption,  the  r-dependent  term  in  Eqs.  (6)  and  (9), 
(1  +  7t)"*,  is  only  s  0.12  for  r  ■  0.1  and  s  0.013  for  r  -  10.  These  numbers  must 
be  compared  with  it  and  with  the  deterministic  shift,  arcsin(a/b),  which  is  0.52.  On 
the  scale  of  our  measurements,  the  predicted  shift  is  just  observable,  but  too  small  to 
measure  quantitatively.  The  quenching  predicted  by  Eq.  (7),  a  a  0.26  for  r  -  O.l  and 
a  a  0.09  for  r  -  lO,  is  also  clearly  evident  and  in  rough  agreement  with  the  widths 
of  the  Pq  curves  shown.  Measurements  for  larger  D  (a  I)  and  larger  a  (-*•  I,  where 
free  running  commences  in  the  deterministic  case)  show  exaggerations  of  all  these 
effects,  however,  the  linearized  theory  cannot  be  expected  to  apply  for  these  large 
parameter  values,  and  so  we  do  not  show  those  data. 

Also  evident  in  Figure  I  is  the  hint  of  a  third  maximum  in  Pq  near  0-0.  This 
is  the  location  of  a  zero  in  the  deterministic  force.  Furthermore  this  zero  is  a  stable 
fixed  point.  In  the  case  of  pure  white,  multiplicative  noise,  the  diffusion  approaches 
zero,  so  that  Po  -*  6(0  -  0)  at  this  point.  This  maximum  is  easily  observed  using  the 
simulator  for  small  t  and  larger  D  values.  The  relative  amplitude  of  this  maximum  is 
observed  to  decrease  with  increasing  values  of  t,  as  the  length  scale  of  the  noise  cor¬ 
related  motions  grows  relative  to  the  size  of  the  deterministic  basin  of  attraction 
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around  0-0. 

In  summary,  we  have  estimated  the  magnitudes  of  the  noise  quenching  and  the 
color  induced  shifts  between  the  center-of-gravity  and  the  maximum  of  a  peak  in  the 
density  using  an  approximate,  linear  analysis  which  begins  with  the  classical  Lan- 
gevin  equation  for  the  CEL.  In  addition,  we  have  qualitatively  observed  in  an  elec¬ 
tronic  analog  all  the  features  forecast  by  the  calculation. 
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encouragement  and  support.  This  work  was  supported  in  part  by  the  U.  S.  Office  of 
Naval  Research,  grant  number  N00014-88-K-0084  and  by  NATO  grant  number 
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where,  for  simplicity,  we  have  assumed  white  noise.  Obviously,  this 
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FIGURE  CAPTION 


Fig.  1.  Three  probability  densities  measured  for  D  =  0.25,  a  -  0.5  and  b  =  I.  The 
values  of  t  are  shown  on  the  curves.  The  vertical  scale  is  in  arbitrary  units,  how¬ 
ever  the  curve  for  r  -  10  has  been  -educed  in  amplitude  as  shown. 


